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MIKHAIL FELDMAN 

Abstract. We consider two variational evolution problems related to Monge- 
Kantorovich mass transfer. These problems provide models for collapsing sand- 
piles and for compression molding. We prove the following connection between 
these problems and nonlocal geometric curvature motion: The distance functions 
to surfaces moving according to certain nonlocal geometric laws arc solutions of the 
variational evolution problems. Thus we do the first step of the proof of heuristics 
developed in earlier works. The main techniques we use are differential equations 
methods in the Monge-Kantorovich theory. 



In this paper we study two models involving limits as p — > oo of solutions of p- 
Laplacian evolution problems. One is a model of collapsing sandpiles proposed by L. 
C. Evans, R. F. Gariepy and the author |[EFG97|| . Another is a model of compression 
molding proposed by G. Aronsson and L. C. Evans \~KE\. In both models the limits 



were characterized as solutions of variational evolution problems related to Monge- 
Kantorovich mass transfer. Solutions that have the form of the distance function to 
a moving boundary arise naturally in the both models. The equations of the motion 
of the boundary for both models were derived heuristically in |[FFG97|| and |[AE|| . 
According to the equations the outer normal velocity of the boundary at a point 
depends on both the local geometry (curvatures) and the nonlocal geometry of the 
boundary. Thus, the motion of the boundary is a nonlocal geometric motion. 

In this paper we prove rigorously the connection between geometric and variational 
evolution problems. Namely, assuming that a moving surface satisfying the geometric 
equation is given, we prove that the distance function defines a solution of the cor- 
responding variational evolution problem. The main assumption is that the surface 
remains convex (or, more generally, semiconvex) during the evolution. For the col- 
lapsing sandpiles model we also prove the corresponding result for the solutions that 
have the form of the maximum of several distance functions (such solutions represent 
several sand cones interacting in the process of collapse). In the proofs we utilize 
the connection between the models and the Monge-Kantorovich mass transfer. This 
allows us to use the differential equations methods in the Monge-Kantorovich theory, 
which have been developed recently by L. C. Evans and W. Gangbo ||EGan|| . 



The examples suggest that solutions of the geometric equations derived in |[EFG97|| 
can develop singularities, even if the initial data is smooth. Thus we do not assume 
below that the moving surface is smooth. This however makes the technique more 
involved. 
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In the forthcoming paper |Fej] we construct convexity preserving viscosity solutions 
of the geometric equations given convex initial data, thus finishing the proof of the 
heuristics developed in |[EFG97]| and \KE\. 

This paper is organized as following. In Sections ^| - [j] we work on the collaps- 
ing sandpiles model. In Section |8] we work on the compression molding model. In 
Appendix |A| we estimate the local Lipschitz constant of the gradient of the distance 
function of a set. 



1. Notation 

Let fi G R™ be a bounded open set. Let T = <9fi. We will use the following two 
distance functions of the set fi. The first, interior distance function, is 

* / \ C distfx, <9fi), xefl, 

The second, signed distance function, is 

m\ js /™\ _ / dist(x,<9fi), iGd, 

lZj ~ \ -dist(x,afi), x £ fi. 

Let x G -R™. We denote Mdn(x) the set of all points of <9fi nearest to x, i.e., 
Nan(x) = {y G <9fi, |x - ?/| = inf |x - z\}. 

Note that Naci( x ) ma J have more then one point. 

Let y G <9fi. Then <9fi is differentiable at y if there exists 5 > such that in a 
suitable coordinate system (xi, ...,x n ) in R n we have 

(3) y = 0, finB { (0) = {x n > ^{x')}nB s (0) where V(x') = o{\x'\). 

Here x' = (x%, x n _i), and Bs(0) = Bg(0) is the ball in R n with radius S and center 
at 0. <9fi is twice differentiable at y if in a suitable coordinate system in R n-1 the 
function \1/ has the expression 

n-l 

(4) y(x') = Y,K l x 2 l +o(\x'\ 2 ). 

i=i 

The numbers ki, K n -i are the principal curvatures of <9fi at 

Let {fit}, t G R+, be a family of open bounded sets in R n . We call a family of 
sets {fit} Lipschitz (with respect to t) if there exists a constant M > such that for 
any t%, i 2 the set (fi tl \ fit 2 ) U (fi t2 \ fitj is subset of both M | t\ — i 2 | neighborhood 
(in R n ) of <9fi tl and M | ti — t 2 I neighborhood of <9fi t2 . 

Let 

£ = U t (fit x {t}) c R n x 

(5) r = u t (on t x {*}) c R n x r^, 

r t = dtt t c R n . 
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Let (y, t) G R n x R 1 ,, be such that y G T t and the function (x,r) — > g?q (x) is 
differentiable at (y,t). We define the outer normal velocity of T t at y as 

(6) V(y) = ^d s nt (y). 

Let (y,to) G T. The surface F is (2,1) differentiable at (y,to) if in a suitable 
coordinate system (xi, ...,x n ) in R n we have?/ = 0, SnBf^O) = {x n > t)} n 
S^ +1 (0) for small 5 > and 

n-l 

*(a/, *) = v(t -t ) KiX 2 { + o(\ x' - y' | 2 + | t - t I)- 

i=l 

£ n is n-dimensional Lebesgue measure. 
7i n is n-dimensional Hausdorff measure. 



2. The collapsing sandpiles model 
In the paper [EFG97| the p-Laplacian evolution problems 



(7) 



d t u p - div(| Du p \ p - 2 Dup) = 0, in R n x (0, oo) 
Up = g, on R" x {0}, 



are considered in the "infinitely fast /slow diffusion limit" p — > oo. Du p denotes the 
gradient of u p with respect to the spatial variables x%, ...,x n . The operator 



A p {u) = divfl Du p \ p ~ 2 Dup) 

is p-Laplacian. The initial data g is nonnegative, Lipschitz continuous function with 
compact support. We assume that 

(8) Lipfo] = L>1 

where Lip[g] is the Lipschitz constant of g. 

The initial data g satisfying @ is unstable in the following sense. It is shown in 
|[EFG97|| that u p u uniformly on R n x (0, oo), and that 



(9) \Du\<l a.e. in R n x(0,oo). 

This limit u does not depend on t for t > 0, i.e., u = u(x). By (H) and @ we get 
u 7^ g. Thus the limit solution is discontinuous at t — 0. The transformation g — > u 
takes place at t — 0. 

The described problem can be interpreted as a crude model of a collapsing sand- 
pile. The function u(x, t) is the height of the pile at the location x, time t. The main 
physical assumption is that a sandpile is stable if its slope does not exceed 1, i.e., if 
Lip[tt(-,i)] < 1. The condition (Rl) implies that the initial profile g is unstable. Ac- 
cording to the model above, the initially unstable sandpile collapses instantaneously 
to a stable one, with the height function u. 

The transformation g — > u was studied in ||EFG97|| by relating it to the evolu- 
tion problem governed by Monge-Kantorovich mass transfer. Introduce the convex 
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functional 1^ : L 2 (R") — > R 1 U {00} defined by 

if v G L 2 (R n ), I |< 1 a.e. 
+00 otherwise. 

Then u(x) = w(x, 1), where w is the unique solution of the evolution problem 



(10) 



j — <9 4 u; G for a.e. t G 1], 

w = $=9, at t = 1, 



with L = Lip[g] > 1. Note that Lip[g] = 1, and thus the evolution problem ( |i0|) is 



well-posed. The connection between this evolution problem and Monge-Kantorovich 



mass transfer is explained in |EFG97]. 



In the static case the limits as p — ► 00 of the boundary value problems for p- 
Laplacian type equations were considered in the papers ||BDM89|| , [|Jan93|| and refer- 



ences therein. In many cases the limits have the form of the distance function to the 
boundary of the domain. 

A limit evolution problem that represents a stable case of the problem (|7|) was 
considered in the paper |[AEW96fl as a model of a growing stable sandpile. Let u(x, t) 



be, as before, the height of the pile, and let f(x,t) > be the sand source function 
which describes the rate at which sand is coming at location x at time t. The dynamics 
is defined as following. Consider the nonhomogeneous evolution problem 



(11) 



d t u p - divfl Du p \ p ~ 2 Dup) = f, in R" x (0, 00), 
Up = 0, on R" x {0}, 



and let p — ► 00. Then for any T > 

Up —>■ u uniformly on R™ x [0, T] 
and this limit u satisfies the equation 
(12) / - d t u e dl^u], for a.e. t > 

with zero initial condition. In particular the case was considered when sand sources 
are concentrated in m points, i.e., if 

m 

/(*,*) = X^WM*)' 

k=l 

where Sd k (-) is the Dirac mass at the point dk G R™. The solution has the explicit 
form 

u = max(0, 24 (i)— | x — d\ \, z m (t)— \ x — d m |), 

where the height functions {^(t)}^ satisfy a system of ODE. Thus the solution has 
the form of growing and interacting cones centered at the points 

By analogy with the static problems and the model of growing sandpiles |AEVV96[. 



one can expect to find solutions of fllHD that have the form of distance function to 
the moving boundary or superposition of such, i.e., growing and interacting cones. 
However, unlike the situation considered in [ AKW9(jfl , interacting spherical cones do 



not maintain their shape under the evolution defined by fllPp. A more suitable class 
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of solutions is the following class of more general "cones" . A single cone solution has 
the form 

(13) w(x,t) = da t (x), 

where {fit} is a suitable expanding family of sets. A multiple interacting cones 
solution has the form (for simplicity we consider the case of two cones) 

(14) w(x,t) = max(d u i(x),d U 2(x)), 

where {fi* }, {fit } are suitable expanding families of sets. Thus in order to find such 
solutions one has to find the appropriate families of sets. 

The following equations of motion of the boundaries of the sets in ( ]T3|) and ([14]) are 
derived heuristically in |[EFG97|| . Assume that T defined by (|5|) is a smooth surface 
in R™ x R 1 . Let y G T t and let k%, k„_i be the principal curvatures of T t at y. 
Denote n = (ki, n n -i)- Let 7(2/) = 7n t (?/) be the radius of the largest ball touching 
T t at y from within fit. Let V(y) be the outer normal velocity of T t at the point y 
defined by (|6]). Then in the case of a single cone solution (|T3| ) the equation for the 
surface T t is: 

V{y) = ^F( K (y) n (y)) for y G T t , 



(15) where 

F(k,j) 



Jo s nr=i (j - K i- s ) ds 



To write equations for the motion of the surfaces = dQ'l, k = 1, 2 in the case 
of the interacting cones solution (|I4"D, we extend the function 7(-) = ■jn(-) of a set 
Q G R™ to a function R™ — > R 1 by setting 7 = on R™\f2, and 7(x) = 7(7/) for 
x E Q such that y G Mdo,{x). Then the equations are: 

V{y) = \F(K(y), ln ,(y), ln}nn 2(y)) for y G T t fc , k = 1, 2, 



rn-l 



/7 s nr=i (! -K%s)ds 



(16) where 

^ ' 7 ' ^ //nr^i 1 ! 1 - 

The equations ( ]T5|) and (p!6|) have the form of nonlocal geometric curvature motion. 
Geometric nonlocality is caused by the functions 7(-)- 

In this paper we prove that if the moving surfaces satisfy geometric equations (|l"5"D 
or ([TBI) an d some regularity conditions, then the functions ( |i~3]) or ( |i~4D respectively 
are solutions of the evolution problem (0). Namely, we prove the following. 

Theorem 1. Let {fit} be a locally Lipschitz continuous family of open convex bounded 
sets, where t G [a,b], b > a > 0. Suppose that the equation f |73[ ) is satisfied at every 
point of (2, 1) differentiability of the surface T defined by (fjj. Then the function (|73[) 
is a solution of the evolution equation 

w 

(17) d t w G dloolw] for a.e. t G [a, b}. 

t 
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Theorem 2. Let {Qj }, {^} be two Lipschitz continuous families of open bounded 
convex sets, where t G [a, b], b > a > 0. Suppose that the equation ftT(\ ) is satisfied 
at every point of (2, 1) differentiability of the surfaces T k = Uj(rf x W)> ^ = 1>2- 
Then the function (\m ) is a solution of the evolution equation (Tip. 



The convexity condition in Theorems [l], ^| can be relaxed to the following "semi- 
convexity" condition. 

Condition 3. (lower curvature bound condition with radius r ). An open 
set Q C R n satisfies this condition if 7 R n\^(x) > 2r for any x G R n \ fL 

It is easy to check that if a set is either convex or has C ,x boundary, then the set 
satisfies the Condition [3[ 

Remark 4. Let the set Q satisfy Condition [|. Let dQ be twice differentiable at a 
point y, i.e., and ^) are satisfied in an appropriate coordinate system on R n . 
Then the principal curvatures of dQ at y satisfy 

Ki > --^ for i = 1, ...,n- 1. 

2r 

This follows from the fact that the graph of the function ^ defined by ([|) and (0j lies 
above the ball -8^.(0? — 2r) by Condition^. 

The following two theorems state that convexity can be replaced by the weaker 
Condition ^ in Theorems [I] and 0. 

Theorem 5. Let {fit} be a locally Lipschitz continuous family of open bounded sets, 
where t G [a, b], b > a > 0. Let for every t G [a, b] the set Qt satisfy Condition [| 
with radius ro > 0. Suppose that the equation ftldj) is satisfied at every point of (2, 1) 
differentiability of the surface T defined by (j^j. Then the function ftH\ ) is a solution 
of the evolution equation (jij 



Theorem 6. Let {^f} be two Lipschitz continuous families of open bounded 

sets, where t G [a, b], b > a > 0. Let for every t G [a, b] the sets flj, Qf satisfy 
Condition [| with radius r > 0. Suppose that the equation is satisfied at every 
point of (2, 1) differentiability of the surfaces T k = U 4 (r^ x {t}), k = 1, 2. Then the 
function ( f/^| j is a solution of the evolution equation ftF\). 

Theorems [I] and ^ follow from Theorems |5| and |6] since convex sets satisfy Condition 
[3] with any radius. 

The proof of Theorems [5] and |B| utilizes the following relation of the equation ( |TTD 
and the Monge-Kantorovich mass transfer problem. Fix t. Equation (|T7|) implies 
that that the function w(-,t) is the Monge's potential for the optimal transfer of the 
measure with density j(-,t) into the measure with density d t w(-,t). Then, according 



to [ EGan|] , there exists a measurable function a(x), the mass transport density, that 
satisfies the following properties: 

a > 0, supp(a) C {x \ \ Dw \= 1}, 

(18) 

— div(a/>u;) = - — d t w, 
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where the last equation is understood in the weak sense and t is fixed. Conversely, if 
there exists a function a(x) satisfying (|I8D , then the function w is a solution of ( [TTD 
for the given t, see Section ^| below. 

The idea of the proof of Theorem [| is as follows. Let a family of sets {^t} satisfy 
fll3|) . An open bounded set can be represented as the union of distance rays, i.e., 
maximal intervals that start at the boundary, on which distance to the boundary 
is the linear function with slope 1. The Condition |3| implies that the collection of 
distance rays possesses "nice" measure-theoretic properties. Fix t and consider a 
distance ray R y of Vt t starting at a smooth point y G dQ t • Introduce the coordinate 
s on R y , the distance to the boundary. The equation (|TJ) can be formally rewritten 
as: 

w 

(19) Da- Dw + aAw d t w = 0. 

t 

Since w(-,t) = dist(-,T 4 ), we have on R y 

, n— 1 

Da ■ Dw = — , w = s, —Aw = > , o t w = V(y), 

as 1 — 

i=l 

where K\, K n _i are the principal curvatures of T t at y, and V(y) is the outer normal 
velocity of T t at y. Thus the equation flT9"|) can be formally rewritten on the ray R y 
as the following ODE 

n-1 

(20) a'(s) - a(s) £ - J + V(y) = 0, 

This ODE has a solution a(s) with zero boundary conditions at both ends of the ray 
R y if V(y) satisfies (|15|). Define a(x,t) by ODE (|20| ) with zero boundary conditions 



on each distance ray that starts at a smooth point of the boundary. We prove that 
a(-,t) is a measurable function that satisfies (|I8| ). In order to do this we examine the 
properties of distance function and distance rays, and use a nonsmooth (Lipschitz) 
change of coordinates on the set Qt with n — 1 coordinates along the boundary and 
the n-th coordinate along distance rays. 

We prove Theorems | and |6] in Sections |3] - 0. In Section |3| we examine properties 
of sets satisfying Condition |3|, and in particular we describe the change of coordinates 
mentioned above. In Section ^ we define the mass transport density function and 
examine some properties of this function. In Section |5] we prove that the mass balance 
equation (|T5| ) is satisfied. In Section || we conclude the proof of Theorem [5[ Finally, 
in Section [I] we sketch the proof of Theorem []. 

3. Properties of distance function and ridge sets of the sets 
satisfying the lower curvature bound condition 

Let Q be an open set. When there is no possible confusion, we write d(x) for dn(x) 
defined by (|). _ 

Let x G Q. Denote by R x the longest line segment through x in Q along which 
dfi(-) is a linear function with slope 1. We call R x the distance ray of x. Note that a 
point x can have more than one distance ray, in such case denote any of them by R x . 
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If x G Q then all such rays R x have equal length. If x G dQ then rays R x in the case 
of general open set Q can have different length. However if Q satisfies Condition |3| 
then there exists at most one R x for x G dQ, see Proposition || below. 

For any xsfi one endpoint of R x lies on dQ and belongs to the set J\fdsi(x). Call 
this endpoint the lower end of R x . Call another endpoint of R x the upper end of R x 
(the names "lower" and "upper" correspond to the positions of the endpoints of R x 
on the graph of y = dn(x)). Let y, z be upper and lower ends of the ray R, and let 
x = Xy + (1 — X)z where A G (0, 1). Then we say that the point x lies in the relative 
interior of the ray R. 

Remark 7. If x G Q lies in the relative interior of a distance ray R x , then the ray 
R x is the unique distance ray that contains x, and the set Afdu( x ) consists of one 



point, see ||EH87 . 



Remark 8. It follows from Remark that sets satisfying Condition [| are sets with 
positive reach defined by H. Federer |Ked59 ]. More precisely, if a set Q satisfies 



Condition [| with radius r then reach(Q)=2r. 

The function 7n(-) in the equations QT5|) , (]Tjj) is called ridge function and can be 
expressed as following: 

x g r™ \ n, 

(21) 7 fi(» = { | R x | x G Q, 

x G dQ. 




The function 7n(-) is uppersemicontinuous in R n , cf. ||tCH87|| , Proposition 3.2. 
Define the set 

TZ = { x G Q | d n (x) = 7n(x) }. 

The set TZ is called ridge set of Q. Note that if x G TZ then either x is the upper end 
of some distance ray of Q, or otherwise x lies at dQ and there are no points y G Q 
such that x G J\fan(y). 

Proposition 9. Let Q C R n be a bounded open set satisfying Condition [| with radius 
r > 0. Let Q r be the r -neighborhood of Q. Then 

(a) 7Rn\n r ( a; ) — r f or a M x £ \ «-e., i/ie sei fi r satisfies Condition^ with 
radius ~. 

(b) For any y r G dQ r there exists a unique y G dQ such that \y r — y\ — r, and 
\y r — x\ > r for all x G fi\{y}. For any y G dQ there exists y r G dQ r such that 
\y — y r \ = r. If dQ is differentiate at y then y — y r is orthogonal to dQ at y. If there 
exists a ball inside Q touching dQ at y, then dQ is differ entiable at y. In particular, 
dQ r is everywhere differentiable. 

(c) For any x G Q we have 

(22) dist(x, dQ r ) = dist(x, dQ) + r. 

If x G Q r \ Q then dist(x, dQ r ) < r and x G' TZ r , where TZ r is ridge set of Q r . 

(d) Sets Q and Q r have the same ridge set, i.e., TZ = TZ r . 
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Proof, (a) Let x G R" \ fl r - Then dist(x, fl) = R > r. Let y G dfl be a point such 
that \ x — y \= R. Let y r be such point of the interval connecting x and y that 
\ y — y r \= r. Then |/ r G fi r . 

Consider first the case R > 2r. Then dist(x, f2 r ) > r. Since if this is not true then 
there exists such point z\ G fl r that | x — Z\ \< r, and there exists z 2 G fl such that 
I z i — z 2 |< r - Now we get i? <| a; — z 2 \<\ x — Z\ | + | z\ — z 2 |< 2r, a contradiction. 
So 

7r™\Tv( :c ) > dist(x, fi r ) > r. 

Consider now the case r < R < 2r. Then by Condition |3| we have 'Jn n \n( x ) — ^ r 
and so R = dist(x, fl) < 7 R ™\n(£)- Thus y is the unique point on dfl such that 
| x — y |= dist(a;, Let z = y + 2r^^f- Then from the condition |3| and |EH87 



Lemma 3.4, it follows that for any x of the form tx + (1 — t)z where < t < 1 the 
point y is the unique point on dfl such that dist(x, dfl) =\ x — y \. Then for each such 
point x, the point y r is the unique point on dfl r such that dist(x, dfl r ) —\ x — y r |. 
For suppose on the contrary that there exists z\ G dQ r such that z\ ^ y r and 
| x — z\ |<| x — y r \. Then there exists z 2 G dQ such that | 2:1—2:2 |< r. Consider 
first the case z 2 ^ y. Then 

I 5 — y I = I x — y r I + I y r — y \ = \ x — y r \ +r > \ x — z\ \ + | z\ — z 2 \ > \ x — z 2 | , 

a contradiction with the fact that y is the unique nearest to x point on dQ. Consider 
now the case z 2 = y. Since the ball of radius R and center x lies outside Q and 
\ z \ — y\ < r i h follows that Z\ does not lie on the interval connecting x and y, and so 

\ x — y \<\ x — Z\ I + I Z\ — y \<\ x — y r \ +r =\ x — y |, 

a contradiction. So, y r is the unique nearest to x point on (9fi r for any x on the 
interval connecting z and y r , so 7 R ™\ ^ (x) >| z — y r \ — r. 

(b) Let y r G <9fi r . Then dist(y r ,f2) = r. Thus there exists ?/ G <9f2 such that 
I V ~ Vr |— r - By Condition ^| , 7 R n\Q(y r ) > 2r > dist(y r , (9f2), and thus by ||EH87|| , 
Lemma 3.4, y is the unique point of dfl nearest to y r . 

Let y G dQ. Let yi G R™ \ Q for i — 1,2, and let yi — > y. By condition || 
ln n \n(yi) — an d by uppersemicontinuity of 7 R n\n(-) we get 7 R n\Q(y) > 2r. Thus 
we have shown that every point of dfl has an exterior tangent ball of the radius 2r. 
So if a point y of dfl has an interior tangent ball, then dfl is different iable at y. 

Let y G <9fi. Then, as we have shown, 7 R n^(?/) > 2r. Thus there exists a distance 

ray R y of the set R" \ fl starting at y such that the length of R y is at least 2r. Let 
e be the unit vector in the direction of R y , let y r = y + re. Then NQ^{y r ) = {y} by 
Remark |7|. Since | y — y r \— r we get y r G 9f2 r . 

We have shown that the interval connecting y and y r lies on the distance ray of 
the set R n \ fl starting at y, and on distance ray of the set fl r starting at y r . It 
follows that the vector y — y r is orthogonal to dfl r at y r and also to dfl at y if dfl is 
differentiable at y. 
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(c) Let x G fl°, y G dfl, and \x — y\ — dist(x, dfl). Then, by (b), 

r 

y r -=y + i -Ay -x)e dfi r . 

\y - x\ 

Clearly, 

\ x ~ Vr\ — \x — y\ + r — dn(x) + r. 
Assume that there exists z r G dfl r such that 

\x — z r \ < dn(x) + r. 

Let z be the point of intersection of dfl with the interval connecting z r and x (such 
point exists since x G fl, z r fl). Then |z — z r \ > r and so 

\x — z\ < dn(x). 

This contradicts the fact that z G dfl. Thus y, r G Aan r (^) and thus (|2|) is proved. 
It remains to prove the last assertion of the statement (c). 

Let x G fl r \fl . Let y r G Afan r ( x )- Let R Vr be the ray orthogonal to dfl r at y r and 
let = Ry r fl <9fi. Then by (b) we have the following: x G R Vr , x lies on the interval 
connecting y and y r , and |y — y r \ = r = dist(y, dfl r ). Then it follows from [EHS7|, 
Lemma 3.4 that x ^ lZ r . Also, dist(x, dfl r ) =\ x — y r \< \y — y r \ = r. 

(d) We show first that 1Z C TZ r . If x G fl and B p (x) C fl, then B p+r (x) C f2 r . Let 
x G 7Z\dfl. Let y G A/an(x) and let y r = y + '"i^Efr- Then by (b),(c) we see that 
y r G <9fi r , and |x — y r \ = dist(x, dfl) + r = dist(x, dfl r ). Let x = y + (1 + e)(x — y) 
where e > is small enough so that x G fl. Since x G TZ, it follows that dist(x, <9fi) < 
1^ — 1/1 = (1 + £)dist(x, <9fi). Let y G <9f2 be such point that dist(x, dfl) = \x — y\, let 
llr = y + r ]frf|5 then y r G c?H r and |x — y r \ = dist(x, dfl r ). Then 

dist(x, dfl r ) =\ y — x \ +r < (1 + e)dist(x, dfl) + r = \x — y r \, 

so 7n r (x) = \x — y r \ = dn r (x). Thus x G TZ r . 

Let x G 7£n <9f2, i.e., there does not exist a ball inside fl that touches dfl at x. Let 
B r (y r ) be a ball of radius r and center at such point y r that B r (y r ) fl O = {x}, such 
ball exists by Condition Then dist(y r , <9f2) = r, and by (b) x is the unique point 
of dfl such that \x — y r \ = r. Let x = y r + (1 + e)(x — y r ) for e G (0, |). Given e, 
there are two possibilities: either x G fl or x ^ fi. 

Let x G fi. Then, since x G 7£ fl <9f2, distfx, dfl) < |x — x| = sr. Let y G <9fi be 
such that \x — y\ = dist(x, dfl), let y r = y + r y|Ef| • Then y r G <9f2 r and by (c) we get 
\ x ~ Vr\ = dist(x, dfl) + r = dist(x, dfl r ), and so dist(x, dfl r ) < (e + l)r = \x — y r \, so 

Let now x ^ fl, i.e., x G fi r \ fl since £ < |. Then, by (c), dist(x, dfl r ) < r, but 
\x - y r \ = (e + l)r, so y r £ M 9 n r (x). 

Thus y r ^ A/an r (x) for any small £ > 0. This implies x G 7£ r . 

Thus Kc1l r . __ _ 

Now we show that ft r C ^. By (b), (c) we get ^ r fl {fl r \ fl) = 0. 

It remains to consider x G f2 fl 7^ r . Let first x G fl 7^ r . Let as above 7/ G A/9q(x), 
yr = y + r ]fzf| e A/"dQ r (x), x = y+ (1 + e)(x — y) where e > is small. From x G 72- r 
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we conclude that dist(x, dfl r ) < \x — y r \ = (1 + e)dist(x,dfl) + r. Let y r G dfl r be 
such that \x — y r \ = dist(x, dfl r ) and let y be the unique point of dfl nearest to y r , 
then \y — y r \ = r (uniqueness of y follows from (b)). Then by (c), (b) the points 
x, y, y r lie on one distance ray and dist(x, dfl) — \x — y\, so 

dist(5;, dfl) = \x — y\ = \y r — x\ — \y r — y\ < \x — y r \ — r = 
(1 + e)dist(a;, dfl) — (1 + e)\x — y\ — \x — y\ 

and so x G 71. 

The case of x G dQ fl lZ r is similar to the case above. □ 



Remark 10. Examples of nonconvex polygons on a plane show that without assum- 
ing some condition of the type of Condition^ the statements (b)-(d) of Proposition 
|] are not true in general. 



Proposition 11. Let Q G R n satisfy Condition^ with radius r. Then dQ p G C 1 ' 1 
for any p G (0, jg). 

Proof. The assertion follows from inequality 4.8(8) of ||Fed59|| and Proposition || (d). 
Also, inequality (|96|) proved in Appendix Al can be applied instead of 4.8(8) of 
Fed5l . □ 



Recall the following facts (see e.g. ||GT83|| , Lemma 14.16, 14.17). Let Vt C R n be 
an open bounded set with C 2 boundary dfl. Let Q r be, as above, r-neighborhood of 
Q. Then dQ r G C 2 for small r. If y G dfl and if y r G dQ r is the unique point on dQ r 
such that \y — y r \ — r then the principal coordinate systems of dfl at y and of dfl r 
at y r are parallel. Denote by Ki and K^r, {i = 1, ...n — 1) the principal curvatures of 
dfl at y and of dfl r at y r respectively. We have 

(23) K-ir = for i = l,...,n — 1. 

1 + Kit 

In the next two propositions we show similar facts for sets with nonsmooth boundaries 
satisfying Condition [5] and for Lipschitz families of such sets. 

Proposition 12. Let fl C R" be an open set satisfying Condition^ with radius r . 
Let < r < r . Fix y G dfl. Let y r G dfl r be such point that \y — y r \ = r. Then: 

a) Let dfl be twice differentiable at y, i.e., in a suitable coordinate system in R" 
we have y = 0, fin B s (0) = {x n > ^J(x')} n B s (0) for small 5 > and ^l(x') = 
Y^=i K i x l + cKI^'l 2 )- Then the point y r is unique given the point y, and in the 
coordinate system introduced above y r = (0, — r), fl r fl Bs 1 (y r ) = {x n > ^/ r (x')} fl 
Bs 1 (y r ) for sufficiently small 5\ > 0, where \l/ r (x') = — r + Yui=i i+l-r ^ ~^ °(\ x '\ 2 )> 



i.e., fl r is twice differentiable at y r and (jl3[) holds, 
b 

at y 
a)). 



b) If dfl r is twice differentiable at y r and 7n(y) > then dfl is twice differentiable 
at y and (1^31) holds (i.e., the functions ^ and \l/ r have same expansions at as in 
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Proof, a) Uniqueness of y r and the fact that y r = (0, — r) follow from smoothness of 
dQ at y. So the interval connecting y and y r lies on x n -axis. By Propositions and 



11 the surface dVt T is of class C ' , and the interval connecting y and y r is orthogonal 



to dVt r . So fi r fl Bs 1 (y r ) = {x n > fy r (x')} fl Bs^yr) for some #i > 0, where \l/ r is a 
C 1,1 function. 

Let e > 0. There exists a > such that 



n-1 



|\P(a/) - ^KjX 2 2 | < e\x'\ 2 if < o. 
i=\ 

Denote A £ = {(z' , z n ) G R n | z n > Y^7=i( K i ~ £ ) z i}- Then we get: 

onnB a (o) c A £ \A-e- 

Since 7 R n^(?/) > 2r, the point y r lies in the relative interior of the distance ray R y 
of the set R n \ Q. Then there exists a\ > such that for any point x r G r\B ai (y r ) we 
have Nani^r) C dfifli^O). (Proof: suppose this is false, then there exist o~i > and 
a sequence rrjf — >• y r such that for x k G Mao,{x k ) we ~ v\ > °i- Then, passing 

to a subsequence, we get x k — > x where x G <9f2 and x ^ y. By continuity of distance 
function we get x G Naniyr) and x ^ y. By Condition ^| y r lies in the relative interior 
of the distance ray R y of the set R n \ Q. Thus y r has a unique nearest point on dQ. 
A contradiction.). So 

(24) x G A £ \A^ £ if s G Afdn(xr) where x r G B ai (y r ). 

Let iV e be r- neighborhood of the set A £ . By Remark |] we have > — Let 
e < ^ Then (by ||GT83|| Lemma 14.16, 14.17), <9iV e is the graph of a C°° function 
9s, ^d ^fH-(O) = .^g^^. Let 

N £ , s = {(z',z n ) | *n>% i+ K :~ £ _ £)r -5)zl-r}. 

i=l 1 ' 

If o~2 > is small enough then 

(25) A> ±£ ,_ £ n B a2 (y r ) C iV ±£ fl B U2 (y r ) C iV ±£i£ n B a2 (y r ) . 
Thus if we choose 03 small, then we get by (|24|), fl25|) 

dn r n B a3 (y r ) cN £jE \N^ e , 
which means that for x r = (x' r , ^/ r (x' r )) G dQ r fl B a . A (y r ) we have 

n— 1 n— 1 

f-' 1 + («i - e)r ' 1 + (k, + e)r 

i=i v ' i=i v ' 

Sending e to we conclude the proof. 

b) By Proposition [5] the condition jn(y) > implies that jn r (yr) > r. Now we can 
perform calculation similar to the proof of assertion a). □ 
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We need similar facts for families of sets. Let {fit} be a family of open sets in R n . 
Let (fit) r be r-neighborhood of fif in R n - Define 

(26) E T = U t [(ft t ) P x {*}], r r = u t [d(n t ) r x {t}\. 

Proposition 13. Let {fit} be a locally Lipschitz continuous family of open sets in 
R n . Let for every t 6 H]_ the set fi t satisfy Condition with radius ro > 0. Let 
< r < r . Fix y G dQ to . Let y r <G d(Q to ) r be such point that \y — y r \ = r. Then: 

a) Let T, defined by ffi), be (2,1) differentiate at (y,t ). That is in a suitable 
coordinate system (xi, ...,£„) in R n we have 

y = o, fin b^\o) = {x n > v(x',t)} n £ 5 " +1 (o) 

for small 5 > and 

n-l 

t) = V(t ~t ) + Y^ K i X l +O{\t-t \ + \x'\ 2 ). 
i=l 

TTien t/ie pomi |/ r is unique given the point y, and in the coordinate system introduced 
above y r = (0, —r), 

fi r n B n +\y r ) = {x n > * r (x',t)} n Bl +1 (y r ) 
for sufficiently small Si > 0, where 

n-l 

$ r (x',t) = -r + v(t - t ) + V - Ki xl + o(\t-t \ + \x'\ 2 ). 

i=i 

b) IfT r is (2,1) differentiate at (y r , to) and 7n tQ (y) > then T is (2,1) differentiate 
at y and the functions ip and ip r have same expansions at (0, to) o,s in a). 

Proof. The proof is similar to the one of Proposition [TJ|. In particular, the sets A £ , 
N £ and N e< s are defined as following. 

n-l 

A £ = { (z', z n , t) I z n > v(t - t ) + y^X K i ~ £ ) z i - £ l £ ~~ *o| }• 

1=1 

N £ is defined as following: for each t = t* the set N £ C\ {t = t*} is e-neighborhood in 
R™ of the set {x \ (x,t*) e A £ }. Finally, 

n-l 

N Ei6 = {(z',z n ) I z B >i;ft-to)+y)( 1 K !~ £ , -S)z 2 -r-(e + 8)\t-t \}. 

^ ± + {Ki- e)r 

The rest of the argument does not change. □ 

To derive further properties of the distance function of a set fi satisfying the Con- 
dition we need to integrate over such sets. To do this it is convenient to decompose 
the set by suitable subsets, and define a local coordinate system in each subset. The 
construction is carried out below and can be roughly described as following. The 
subsets are unions of distance rays passing through subsets of the boundary. The 
coordinate systems consist of variables x\, ...,x n -i on the boundary, and the variable 
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x n along the distance rays. However, since dQ is not smooth enough, we have to use 
dQ r . Now we turn to the construction. 

By Proposition 11 the set fl r for small enough r > has C ' boundary. Fix such 
r. We can choose sets U\ } ...,Un, bounded and open in R n , so that dtt r C Uf =1 Z4, 
and for each k = 1, iV in a suitable coordinate system in R n we have 

Q r nu k = { (x',x n ) | x e U k , x n > $ k (x') } nu k , 

where U k — { x' \ (x 1 , x n ) G U k } C R" _1 is an open set, and <& k is a C 1 ' 1 function on 
R nl . 

We can also choose open sets Uq and Un+i C R n such that Uq C VL r , Un+i C 
R n \f2 r , both Uq and U N+ i do not intersect <9f2, and R n = U^ 1 ^- Let k = 
1, N, be a smooth partition of unity on R n related to the sets Uo,Ui, ...,Un+i, i.e., 
J2k=o = 1 on R n and & k G C^(U k ). Then by our choice of U and Un+i we get 
Z k= i ^k = londa r . 

Define maps G fe : U k x R 1 -> R™ (k — 1, JV) by 

(27) G fe (x',x n ) = y + x n Dd s nr (y) 

where y = (x', $fc(x')) G <9fi r , and c?Q r (-) is the signed distance to dfl r . 

For each z G <9fi r the vector Dd^^z) is the inner unit normal to dQ r in z. Since 
<9f2 r is a C 1 ' 1 manifold, -Dc^ r is Lipschitz on dQ r . Then since is C 1 ' 1 on R n_1 , 
the map G k is Lipschitz on bounded subsets of U k x R 1 . In particular, the map 
U k — > dQ r defined by x' — > G k (x', 0) is Lipschitz. 

From that by explicit computation we get the following 

Lemma 14. Let the map G k : U k x R 1 — > R n be defined by Then the Jacobian 

JG k is a locally bounded measurable function. For £ n_1 a.e. point x' G U k the 
surface dfl r is twice differentiable at the point G k (x',0). For every such point x' the 
map G k is differentiable at (x',x n ) for every x n G R 1 . At every such point (x',x n ) 
the Jacobian of G k is given by 

n-1 

(28) J n G k { X ', X n ) = + \D$ k { X ')\ 2 J] [1 - K r>i X n ] , 

i=l 

where the principal curvatures of dQ r at G k (x',0). 

Denote by U k the sets 

(29) U k = U x , et j k { (x',x n ) | x n G [0, 7nr (y)) } C R", 

where 7o r : R" — > R 1 is the ridge function (^1|) of the set fl r . Note that jn r (z) > r 
for all z G fl r by Proposition [5] (c, d). 



The map G k is one-to-one on U k (by definition of 7n r and Lemma 3.4 of |[EH87|| ). 
In addition Q\7Z C ^ k=1 G k (U k ), since every point of Q r \7Z lies inside some distance 
ray. For each k — 1, N define a function ^> k : R™ — > R 1 by extending the function 
ty k from the boundary inside Q r as a constant along distance rays and defining ty k 
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to be the on the ridge of Q r and outside Q r , i.e., 

(30) * (z) = { * fc( ^)' where y G N 9 n r (z), if z G O r \ K 

\ 0, otherwise. 

Note that y is uniquely defuned by z in the first case of ( |30|) . We have 

N 

supp^ fc C G k (U k ); J2^ k (y(z)) = 1 for all z G tt r \1l. 

k=i 

We thus decomposed the set Q r \ 1Z by the sets G k {U k ) and defined a related 
partition of unity ^ k . Note that in general the sets U k and G k (U k ) are neither open 
nor closed, and the functions ^ k are discontinuous. In the next propositions we prove 
that these sets and functions are measurable. 

Proposition 15. Define a function j r>k : R n_1 — > R 1 by 

/o-n % (~r\ - S 7n r (G k (x',0)), if x' G cl(U k ) 

{S) lr > k{X) -\0, ifx'icl{U k ) 

where cl(-) denotes closure of a set. Then 7^ is uppers emicontinuous (use). The 
sets U k are C n measurable. The ridge set of Q r has C n measure 0. 

Proof. It is enough to show that for any sequence {x' m }™ =1 such that x' m G cl(£4), 
x' m — > x' and ^ r ,k(x f m ) — > L we have 7 r ,A ; (x / ) > L. 

The function 7n r is use on R n . Since x G cl(U k ) and G k is Lipschitz on c\(U k ) x 0), 

%k(x') = 1 n r {G k (x',0)) > hm jn r {G k ({x'J). 



m— >oo 



Thus, 7^ is use. 

We have [4 = { ( ' G C/fc, < x n < 7r,fc(x / ) }. Thus C/jt is ^"-measurable 

(note that U k is open in R n_1 ). 

To prove that the ridge has C n measure 0, we note that G^ 1 (!Z)nU k = {(x',x n )\x' G 
U k , x n = 7 rj fc(x')}. Then £ n (G A T 1 (7?.) Pi U k ) — because J r>k {-) is use. Since G k is 
Lipschitz on c\(U k ) x [0,diamf2 r ] and since D k=1 U k = f2 r , we have C n (1Z) = 0. □ 

We will use repeatedly the following fact: 

Lemma 16. Let f : A C R m — > R n (m < n) be a Lipschitz function, A be a C m - 
measurable set, and let g : R m — > R 1 6e a C m -measurable function. Then the function 
g o : R n — > R 1 C n -measurable. 

Proof. We can assume that g = outside A. Let K C R 1 U {+00} be an open set. 
Then W = g~ l (V) n A is £ m -measurable. Then (# o J" 1 )"^ = f{g-\U )) PiA = 
f{W), which is a measurable set since / is Lipschitz and m < n (see [|EGar92| . 
Lemma 2 of 3.3.1). □ 



Proposition 17. Functions ty k and ty k o G k are C n -measurable. 
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Proof. The functions G k are Lipschitz. Then by Lemma [16] it is enough to prove 
that functions ty k ° G k are £ n -measurable (since then ty k = (^ k o G/-) ° G^ 1 ). Let 
(x', x n ) G U k . Then 

V k (G k (x',x n )) = * k (G k (x', 0)) = y k (x',$ k (x')). 

Consider the function Q(x',x n ) = ^ k (x',^ k (x')) . Since ^ k G C°°(R n ) and $ fe G 
Ci,i(R»-i) we see that 9 is a C 1 - 1 function R n -> R 1 . Now since f fc oGj = 9on 



C/fc, ^oGt = outside £4, and by Proposition jl^ the set U k is measurable, the 



assertion follows. □ 

Proposition 18. Let Q be a bounded open set satisfying the Condition [|. Let X C 
<9f2 and 7i n ~ 1 (X) = 0. Lei F 6e the union of distance rays passing through X , i.e., 
Y = U X€X R X . ThenC n (Y)=0. 

Proof. We can assume that X C G k (U k ) D <9f2 for some fc, and then F C G k (U k ). We 
also can assume that 

(32) 7n(z) > for all x G X 

If ([32]) is not true then we can replace X by X n {x \ > 0} and F does not 

change. 

Since j rjk (x') is a u.s.c. function of x', we get C n ({x n = 7^(2/, 0)}) = 0. We 
also know that J n Gfe > on the set U k . The map is one-to-one on U k , and so 
N(G k ,y) = 1 for w G F C G k (U k ), where N(f,z) denotes the multiplicity function, 
which is the number of elements of f~ l (z). Then by area formula we obtain 



(33) C n (Y)= XY (y)dy= X Y(G k (x))J n G k (x)dx : 

JR" Ju k 

where xa(') is the characteristic function of the set A, i.e., xa(-) equals 1 on A and 
outside A. By Lemma III the Jacobian J n G k is a locally bounded measurable 



function. Let x' G U k be such point that the surface dVL r is twice differentiable at the 



point Gk(x',0). By Lemma [14] the Jacobian J n G k is defined by the expression fl28|) 
at (x',i„) for every x n G R 1 . The principal curvatures of dfl r at the point G k (x', 0) 
satisfy 

(34) K ri > --, «r,i7r,fc(z') < 1 

r 

by definition of 7V,fc(-)j an d ^y Condition |3|. We also have 

(35) Tv^a/) < diamf2 r . 

Let m < n be such that — - < K r i, /t r ?n < and fiv jm +i, K r,n-i > 0. Then for 
x n G (r, 7 r! fc(x')) we calculate using (|34"D , (p5|): 



(36) < JnGkjx 7 ^n) ^ rit=l [1 ^r.^ra] ^ 



I + -7r-,fc(^) 

r 



<c, 



where C depends on r, n and diam(O). By Lemma |TJ| the inequalities ( |3"6"D hold for 
a.e. x G f/fc fl {(a/, x n ) | x n > r}. 
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Define a map G k : U k -> dtt C R n by C7 fc (x') = G fc (x', r). Then by (|8]) 
(37) J n _iG k (x') = J n G k (x',r) for a.e. x' E U k . 

So using the fact that the relation Y <Z Q implies XY{G k (x',x n )) = for x n < r, we 
get from ([33|), (0), (0) and (j3"?D using area formula: 

C n (Y) < C [ X y(G k (x)) J n ^G k (x')dx 

= C J n -iG k (x')[ XY(G k (x',x n ))dx n }dx' 

JU k Jr 

= C [ xx(G k (x'))[% k (x')-r]-J n ^G k (x')dx' 
< G xx(G k (x')) J n - 1 G{x')dx 



= C X x(z)N(G k ,z)dH n - l (z)=Cn n ~ l (X) = 0, 
Jan 

where we have used the fact that N(G k , z) = 1 for z G X (which follows from 
and Proposition ^ (b,c) ). □ 

Remark 19. Without assuming a condition of the type of Condition^ the assertion 
of Proposition 18 is not true. For example consider a nonconvex polygon in R 2 and 
take X to be one point, a vertex of a re-enterant corner. 

Remark 20. Note that it follows from the proof above that 

C n {G k \Y))= f J n ^G k (x')i r MX) X Y(G k (x',x n ))dx n ]dx' = 0. 

The following proposition describes the structure of the boundary of a set that 
satisfies Condition [3[ 

Proposition 21. Let Q be a bounded set. Let Q satisfy the condition^ with radius 
r . Then dQ is (Ti 11 " 1 ^ — 1) rectifiable subset ofH n . In addition, dQ is twice 
differentiate H n ~ 1 a.e. on the set 

(38) B = {x E dVl | 7n(a;) > 0}. 

Proof. Let r > be such number that dQ r is C 1,1 . Existence of such r follows from 
Proposition 11. It follows from Proposition |9] b) that the nearest point projection 
mapping P : Q r — > Q is well-defined and onto. Moreover, by inequality 4.8(8) of 
Fed59|| (or by inequality (|95| ) proved in Appendix [A| below) the map P is Lipschitz. 
Since dQ r is C 1 ' 1 , we get that dQ is (7i n ~ 1 ,n — 1) rectifiable. 

By Proposition |9] every point of the set B has a unique nearest point on dfl r . Since 
dfl r is C 1 ' 1 , it follows that dQ r is twice differentiable 7i n ~ l a.e. Since the map P is 
Lipschitz it follows that for 'H n ~ 1 a.e. y G B the surface dfl r is twice differentiable at 
the corresponding point y r . Applying Proposition |T2| b), we conclude the proof. □ 
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Propositions [15|, [18| and [H] imply 

Corollary 22. Let fi C R n satisfy Condition [|. Then C n a.e point of Q lies in the 
relative interior of a distance ray that intersects <9fi at a point at which <9fi is twice 
differentiable. 

We need an analogue of the Proposition [H] for a Lipschitz family {fit} of sets that 
satisfy Condition |3| with a radius r independent of t. Let E, T be defined by (0). 
Let E r , T r be defined by (|26|). Denote 

(39) D(x,t) = d s nt (x), 

where dn t (x) is the signed distance function (|2]). 

Proposition 23. Let {fit} be a locally Lipschitz family of open bounded sets that 
satisfy the Condition [| with radius r independent oft. Then 

a) C n+l a.e. (xo,to) £ E lies in the relative interior of a distance ray of the set 
fi fo that intersects dfl to at such point y G dQ to that Y is (2,1) differentiable at (y,t ); 

b) There exists a constant C depending only on ro such that the function D(x, t) — 
C\x\ 2 on the set E ro is concave in x for every t. 

Proof. First we prove the following fact. 

Lemma 24. Let fi satisfies Condition [| with radius r . Then the distance function 
dn(x) is semiconcave in fi, i.e. there exists a constant C depending only on r such 
that the function d^(x) — C\x\ 2 is concave in fi. 

Proof. Let i6fl. Let w G R/ 1 satisfy 

M < d n (x). 

Let y G Afdn ro (x), where Q ro is r - neighborhood of fl. Then by Proposition ||(c), 

\x-y\ = d n (x)+r , 
\x + w-y\ > d n (x + w)+r , 
\x — w — y\ > dn(x — w) + r . 

Then we have 

C 

da(x + w) — 2c?n(x) + dn(x — w) < \x + w — y\ — 2\x — y\ + \x — w — y\ < — \w\ 2 , 

r 

and so d^ — 2 — \x\ 2 is concave. □ 

Now we prove Proposition [23|. The function D(x,t) defined in (|39|) is Lipschitz 
since {fit} is a Lipschitz family of sets. By Proposition |^ the set (fit) ro satisfies 
Condition ^ with radius y, and 

d (n t ) ro ( x ) = d h t ( x ) + r o for ( x > f ) e Era- 
Then by Lemma the function D(x,t) — C\x\ 2 on the set E ro is concave in x for 
every t. Now by Theorem 1 of Appendix 2 of ||Kry87|| the function D(x,t) is (2,1) 
differentiable at C n+ a.e. point of E ny It follows from Proposition ^ that if D(-) 
is (2,1) differentiable at a point (xo,t ) G E, then T is (2,1) differentiable at (y,to), 
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where y G J\fan t ( x o)- Note that x lies in the relative interior of its distance ray if 
dn to is twice x-differentiable at Xq. The Proposition is proved. □ 



4. Properties of mass transport density. 

Let O be a bounded set that satisfies the condition with radius r . Let t > 0. 
Define a function a = oq^ : R n — > R 1 as following. Let x G dfl be such a point 
that the surface dfl is twice differentiable at x and let Ki, ft n _i be the principal 
curvatures of <9fi at x. Introduce the coordinate s on the distance ray R x by s(z) = 
\z — x\ for z E R x . Then s changes in the interval (0, 7(2?)). Define the function a on 
as the following: 



n— 1 -, /.o n— 1 



1 ' _ 1 /" s _ _ 



i=l * i=l 



/ Q 7(a; ^m=i 1 (i-M^^ e 



/ 7(s) nr=t(i-M^ 



By Corollary ^2] the function a is now defined a.e. in fl. 
Define a = on R n \f2. Now a is defined a.e. in R n . 

Definition 25. The function a(-) = an,t( - ) 2s called mass transport density. 

Define the function V : dQ — > R at the point of twice differentiability of dfl as 
following. 

( 41 ) v ( x) = i^n^'d-^ 

Proposition 26. Lei f2 G R" fre a bounded open set that satisfies Condition [| with 
radius r . Then there exists a constant C depending only on r , n, diamQ such 
that the following is true. Let x G dQ be such point that dfl be twice differentiable 
at x. Then the mass transport density on the distance ray R x , defined by ^J\ ) for 
s G [0,7(0;)], satisfies: 

(a) \a(s)\ < C . Derivative a'(s) exists and is continuous for s G (0,7(2;)). 

(b) a(s) satisfies the equation 

n-l 

^)-«(.)En 7+V(x)=0 

1 fciS b 

1=1 

on (0, 7(2:)), the boundary conditions a(0) = 0, 0(7(2;)) = 0, and a(s) > on 
(0,7(2;)). 

fcj For s G (0,7(2;)) £/ie inequalities hold 

\a{s)\ < C(j(x) -s), |a(s)|<C7s. 
(dj |o'(s)| < C on (0,7(2;)). 

Proof, (a) Existence and continuity of a'(s) on (0,7(2;)) are checked explicitly. The 
bound \a(s)\ < C will follow from (c). 



20 



MIKHAIL FELDMAN 



(b) The equation is checked explicitly. The boundary conditions follow from (c). 
The inequality a(s) > on (0,7(2;)) holds since either «j < 0, or < 7(2;) < — and 
then, defining 



V(s) 



we calculate: 

rn-l, 



dV s ULii 1 ~ M Jo nili^ 1 - Kirfdu - Utii 1 - K i s ) Jo v ULii 1 - w)du 



ds 



[/;nr=i(i-^w] 2 



nr=i 1 ( i - 



10 

n—l 



[/ s nr=i(i-M^)] 2 io fi 



/ IT^ 1 _ K i u )( s ~ v )dv > for s G (0,7(2;)), 



and so 



1 n-l -J 

s ) = -[V r ( 7 (x))-y(s)] JJ— / J^l-ze^tfe/X) for se (0,7(2;)). 

t . =1 1 Jo , =l 



(c)We write 7 for 7(2;): Denote $(s) = /„' ^IiSC 1 -^u)du, V(s) = J* U^ii 1 
Kiv)dv. For s G (0,7) we calculate: 

M-yMs) - $( 7 )*(7)] + MlMl) - $(s)*(7)] ^ _ 



ict(s) 



(7-s)- 



*(7) 

-$(7)^) + $'(s 2 )fr( 7 ) 
*(7) 



n—l 1 
i=l 



n—l 

n 1 - KiS : 

i=l ' 



where si,s 2 G [s, 7]. Note, that V(£) = ^tIt. So 



ta(s) 



n—l 

1=1 



n—l 



n—l n—l 
8=1 



1=1 



(7-s) 



A 7 1 - 1 1 - 

i=i i=i 



n-l 



1 - KiS 2 



(7-s). 



If Ki > then for s* G (s, 7) we have < SKi < s*Ki < 1, and so 

1 — KiS* 

< — < 1. 

1 — KiS 

If Ki < then it follows from the Condition 13 that < C, and so 



< 



1 + 1 


Kj\ 


S* 


1 + 




s 



< 1 + C7 < 1 + Cdiamfi. 



It remains to estimate V(s) for s G [0,7]. 

rn— 1 

J "l 

(42) 



V ( s ) < rs ^t-l: rr~ = 7 < diamO. 



/o" nr=i (! - 
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So, ta(s) < Cdiam fl (7 — s), and, since t > 0, 

(43) < a(s) < C7(7 - s) for s G (0, 7(2;)). 

To estimate a(s) near s = 0, we compute using ([42"D 

|ta(s)| = 



< 



ra— 1 1 p S n—1 

[V(l)-V(s)}H- / H(l- Ki v)du 

i=i 1 KiS - 70 i=i 

n—1 -| „ s n—1 n—1 ^ 

i=i * >/u i=i i=i 



1 — K,S* 



KjS 



where s* G [0, s]. Since ^7 < 1, then for s < ^ we get 

1 — KjS* 



< 



1 — K4S 



< 2 if > 0, 



and since k,- > — - we get 



< - — — < C(diam(Q), r) if «< < 0. 

1 — KjS 

So < C7S, and, since t > 0, we get 

|a(s)| < Cs. 
If s > 2' t nen s — 7 — s > an d from (|43l) 

|a(s)| < C(7-s) < Cs. 



(d) We prove that \a'{s) \ < C for s G (0,7). 



n-l 



a s = as 



^ 1-KiS t 

i=l f 



By ([42"D, [ V(rc) | < j. It remains to prove that 



(44) 



n-l 



a(s 



1 - K 



i=l 



< C 



If > then < ^7 < 1, and since 7 > s we get: 



< 



ft-; 



< 



1 — KjS Kj7 — KjS 7 — s ' 



If Kj < 0, then 



1 — KaS 



1 

< -. 

s 



Now (041) follows from (c). 
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5. Mass balance equation 

Let {fit} be a Lipschitz family of open sets satisfying Condition |] with radius 
r > 0. Let {fit} satisfy equation flip] ) on the time-interval [a, b] in the sense described 
in Theorem |5|. Let the function w(x, t) be defined by (|X3|) . Let a(x, t) be the function 
a n t A x ) fr° m Definition |25. 

The purpose of this section is to show that for a.e. t G [a, b] 

(45) / aDwD(pdz = (— — w t )<pdz 

Jr™ Jr." ^ 

for any G C°°(R n ), where u> t denotes c^u>. Thus we show that the function a(-,t) 
satisfies (|T^) for a.e i. 

It follows from Proposition |23| that for a. e. t the following is true: C n a. e. point 
x G fit lies in the relative interior of the distance ray that intersects <9fit at such point 
y G <9fi t that the surface Y is (2,1) differentiable at (y,t), where T is defined by @. 
Fix such t. For the rest of this section we drop t in the notation, i.e., we write fi, 
a(x), w(x), w t (x) for fi t , a(x,t), w(x,t), w t (x,t). 

We use same notation a(-) for both mass transport density a(x), a function defined 
on R n , and for mass transport density on a ray R x , the function a(s) defined on 
[0, 7 (x)]. 

Denote by T the subset of <9fi that consists of all points at which the surface T is 
(2, l)-differentiable. By choice of t and Proposition |21] we get 

7T _1 (B\f) =0, 

where the set B is defined by (^8]). From definition of B it follows that no distance 
rays of fi have their lower ends in the subset <9fi \ B of the boundary. Denote by 
-R(r) the union of all distance rays that have lower ends in T, i.e., -R(r) = U xe fR x . 
Then it follows from Proposition [T8| and Remark ^ that 

(46) £ n (fi \ R(T)) = 0, C n [G^(n\R(r))} = for k = 1,...,N. 

Let V(-) be the function (^Tj). It follows from that at a. e. point a; of fi 
the function w t (x) is given by the expression \V{y) where y G Ndsi(x). In the next 
proposition we prove that w t is measurable. 

We continue to use local coordinate systems on fi r , described in Section |3| and 
notation introduced there, in particular sets £4, partition of unity ^ and coordinate 
mappings defined by (|29|), (|30D and (|2^) respectively. 



Proposition 27. a) For k = 1, ...,n t/iere exists a bounded C n measurable function 
Vk '■ Uk —> R 1 such that 

(47) v k (x',x n ) = V(G k (x',r)) if G k (x',r) G f, x n G [r, 7 r ,fc(x')]. 

5j T/ie function w t is C n measurable and for any bounded measurable function 
if. Uk^R 1 

(48) ^ / v k tpJ n G k dx = I w t (<p o G k l )dx 



t 



U k n{x n >r} JG k {U k ) 
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Proof, a) dfl r is a C 1 ' 1 manifold, and so the second fundamental form of dQ r is de- 
fined as a differential form with L°° coefficients. At all points where dfl r is twice 
differentiable the principal curvatures of dfl r are the roots of the characteristic poly- 
nomial of the second fundamental form. To write this in coordinates, we recall that 
dQ r DGk{Uk) is the graph of a C 1 ' 1 function x n = defined on the set Uk G R n_1 . 

Then there exist L°° functions [uij{x')]™ such that at the points where -D 2 ^ exists, 
„.. 

dxidxj t 3 

Denote by P x '{s) the characteristic polynomial of the matrix [wy(x')] i.e. P x >{s) = 
det([uij(x')] — s/ n _i). Note that the coefficients of the polynomial P x '( 4 ) are mea- 
surable functions of x' since they are compositions of polynomials and measurable 
functions Uij(x'). 

The eigenvalues of the matrix [uij(x')]™j =1 are principal curvatures of dfl r at the 
points of twice differentiability of Let x' be such a point, let K rj ±, K r)n _i be 
principal curvatures of dfl r in Gk(x',0), and let m be such number that ^ for 
i = 1, m and K r>i — for i — m+ 1, n— 1. Then P x '(s) = [T^Ti ( K r,i — s )) an d so 

n—1 m m _ 

( _i)-v-ip x ,(_) = n (i - «r,is) = n^ri( — *)■ 

i=l i=l «=1 ' 

For any x' G £4 the point G^x', r) belongs to dQ and 0) — Gk{x', r)| = r. If 

the surface <9f2 is twice differentiable at the point y = Gk{x',r), then, by Proposition 
O the surface <9f2 r is twice differentiable at the point y r = Gk(x', 0), and, denoting 
by {/tj}™ =1 and {/t rj i}™ =1 the principal curvatures of <9f2 at y and of <9£l r at y r , we get 

so = h r (if 7^ 0). 



Then, assuming as before that «i, n m ^ and K m+1 , K n _i = 0, we get if m ^ 



n-1 „. „. „. 

na-M = n^nt- s )=rhn 

i=l i=l i=l i=l i=l 



(s + r) 



t— rm -| 
1 lj = l K r,i s + r 

If m = (i.e., all ^ = 0), then - «»*) = (~l) n ~\s + r^P^^) = 1. 

Now we see that at all x' as above, the function V defined in ([II]) can be expressed 
at the point Gk(x', r) G dVL as 

since, by Proposition |9|, 7n r ,(z) = 7(2)01 + r for z G O. 

Define a function 77 = r](x',x n ) on t/jj. x R 1 by the right-hand side of (|I9D. This 
function 77 is ^"-measurable: r\ does not depend on x n , the coefficients of the polyno- 
mial P x '{ 4 ) are measurable functions of x' and by Proposition |TS|, the functions 7^ 
are measurable functions of x' . 
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Then it follows from Proposition [15| that the function 

v k (x',x n ) = rj(x' , x n )xu k W , x n ) 

satisfies all properties asserted in (a) (note that the right-hand side of (|47|) is bounded 
by diam Q). 

By Lemma [16] the function v k ° G^ 1 is £ n -measurable. By (|46|) we get 

\v k o G" 1 = w t at C n a.e. point of Gk{Uk H { x n > r }), 

(50) 

and «; t = on GkiUk fl { x n < r }) 

and so w t is measurable. 

Since map Gk is one-to-one and Lipschitz on Uk, and since Uk is measurable, we 
get by Theorem 3.2.5 of |[Fed69|| and Lemma [16] that for every bounded measurable 
function if 



(51) / v k (pJ n Gkdz= (v k oG k 1 )(x)((poG k 1 )(x)dx 

JU k JG k (U k ) 

Now (M) follows from (0), (M). □ 



Corollary 28. The function j — w t is measurable, and for each ip G C°°(R n ) 



(52) / (y -w t )tpdz = - y~] / {woGk-Vk)((poGk)(^k°Gk)J n Gkdx. 

7r« t t fc=1 JG-^RiT)) 



Proof. We have 



AT 
k=l 



Since supp C Gfc(t4), the integrand on the right-hand side of ( |52"D is measurable 
by Propositions [?7| and [Tj| It is also clear that the integrand on the right-hand side 
is bounded. 

Then from Proposition |27| (b) and area formula, we get 



R" 



1 N f 

)ipdz = - V] / (woG t -tj t )(p GfeX^fc o G k )J n G k dx. 
f k=i jRn 



N 

,W 



This formula is true for each ip £ C°°, and, by approximation, for every bounded 
measurable <p. Now (06|) implies (B2p. □ 



Now we prove that the mass transport density a(-) is measurable. The proof is 
similar to the proof of Proposition ETl 



Proposition 29. The function a(-) is C n -measurable function on R™. The function 

k 



ooGj, 1 is C n -measurable function on Uk- 
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Proof. For (V, s) such that 

G k (x',r) G f , s6 [r,7 rjfe (x',r)], 

the function a o G k (x', s) is given by expression (|40|) , where the curvatures are com- 
puted at the point G k (x',r). Calculations similar to those in Proposition ^7] show 
that this expression can be rewritten using the characteristic polynomial P x >{ m ) of the 
matrix D 2 Q k { following: 



a o Gu(x', s) 



t ( s + r )n-ip x ,(-L-) 



X 



(53) 



r 

Jo 



[v + r) n 



u+r 



By (M), the function a o CT ] 



)dv ^- r (u + r)^P xl (^)du 



• k y j is defined by the above formula at a.e. point of U k , 
and is £ n -measurable function of variables x', s on U k since it is a rational function of 
variables s and jk,r{x') with measurable coefficients and jk,r(%') is C n ~ l measurable 



function of x' . By Lemma [IB, the function a(-) is measurable. □ 
Proof of ( fl5|) . We will transform the right-hand side of (|52|) . From (|30|) we get 

V k (G k tf, s)) = y k (G k (x', r)), sG (r, % k (x')). 

dQ be the map defined by 

G k (x') = G k (x',r). 

Then for k — 1, N, using definitions of u>, v k , ^ k we compute 
1 



(54) 

Let G k : C/fc 



t 



(w oG k - v k )((p o G k )(^ k o G k )J n G k dx 



(55) 



G fe 1 ( r ) 



* fe (G fc (x',r))^ / -[s-r 



v k (x, s)]<p(G k (x' ', s))J n G k (x', s)ds \ dx 



Let a;' G R n 1 be such that G k (x',r) G T. Then <9f2 is twice differentiable at 

G k (x', r), and so 

^f fc (x',s) = ^(L7 fc (x',r)) for s G [r, % k (x')}. 

By Proposition |12], <9f2 r is twice differentiable at G k (x', 0). Let K n _i be principal 

curvatures of <9fi at G k (x',r), and fc rj i, K rjn _i be principal curvatures of dfl r at 
Gfc^O). Then K ri j = ^ k . Let ?/ = G k (x',r). Then the mass transport density on 
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the ray R y satisfies ODE from Proposition pB|(b). Thus we get for s G [r, •y r ,k{x')] 
[r, 7n(l/) +r]: 



dais — r) 



n-l 



,s-r- v k 



+ a(s 



(is ' 1 — Ki(s — r) 

da(s — r) 



'/- 



8=1 



where a(s) is mass transport density on the ray R y . 

By Proposition ^ exists and is continuous and bounded on the interval (0, Jr tk (x')- 
r), and 

We also see by (|28|) that 

(56) ^(x', S ) = J n G fc (a;', S )^ r ^ 



a(0) = a(7r lfc (a/) - r) = 0. 



n-l 



<9s 



Then we calculate, integrating by parts: 



~[s - r - v k (x', s)]<p(G k (x', s))J n G k (x', s)ds 



da(s — r) 



ds 



a[s — r 



q s JnG k {x i s) 

J n G k (x', s) 



(p(G k (x', s))J n G k (x\ s)ds 



a ( s ~ r)-^-(f(G k (x', s))J n G k {x', s)ds 

OS 



This allows to derive from (52[), (|55|) and (p4|) the following equality: 



(57) 



— - w t )(fdz 

EL ^ 



fe=1 JG^<?) Jr 



d 



^k{G k {x', s))a(G k (x', s))—tp(G k (x', s)J n G k (x', s)dsdx' 



Let z = G k (x',s). Using the fact that distance function is twice differentiate on 
R(f), we get: 

^-(y{G k {x',s)) = (D<p)(G k (x',s))-j-(G k (x',s)) 

(58) = Dip{z) ■ Ddn(z). 

We insert the right-hand side of ( |58| ) into the right-hand side of ( |5TD and use area 
formula and ( |6| ) to change variables from (x',s) to z in the integral on the right- 
hand side of fl5"T|). Then using that ^ fc ^ = 1 on Q°\TZ and that a = on 7?. and in 
R" \ fl, we get 
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6. Proof of Theorem |5] 

To conclude the proof of Theorem ^ we need to demonstrate the following. Let, as 
above, w(x) = dn(x), and let v G W /1 '°°(R ra ) satisfy \Dv\ < 1 a.e. Then 



(59) / (^-w t )w> / (y-tw t )v. 

We show that this follows from fl45|) . By approximation, fl45|) is true for Lipschitz 
Let ip = w — v. Then by (^) that it is enough to prove that 

/ aDw(Dw - Dv) > 0. 
Jr.™ 

But, since a > in R n and a = outside Q° and on 7Z, and |£)w| = 1 on Q\1Z, 

/ aDwDw = I a > aDwDv 

JR n JR" Jr» 

since \DwDv\ < ^\Dw\ 2 \Dv\ 2 < 1. 
Theorem |5| is proved. 

7. Proof of Theorem |6] 

The proof of Theorem |6| is similar to the proof Theorem |5|. We give a formal 
calculation. Each step can be justified the same way as in the proof of Theorem |5|. 
Fix t G [a, b}. Define 

V l = v x (t) = {xe n\ | d n i(x) > d^(x)}, 



2 



V 2 = v 2 (t) = {xen 2 t \ d n 2(x) > d a i(x)}. 

Let r > 0, let Ql = fij r and Q 2 . = Vt 2 tr be r-neighborhoods of Q 1 = ft] and Q 2 = Q t , 
and let r be so small that dfll and dQ 2 are C 1 ' 1 manifolds. 

Let 1=1,2. Let {C^}^ and {U l k } k l =l be the coordinate neighborhoods associated 
the set tt r same way as {Uk}, {Uk} with Q r in the proof of Theorem ^ Let G l k : 
U k — > R™ be the corresponding coordinate mappings (0) and let be the 

corresponding partitions of unity (p0|). Let 7' fe be the functions defined as in fl3T|) 



from the functions 7^; using the coordinate mappings G k for k = l,...,Ni. Define 



r k 



the following functions 5 rk for k — 1, ...,iVj. Let £(•) = 7Qinc 2 (")- Then 

# = / 5(Gt(x',r)) +r, if x' G cl(^) and Gj^r) G d^ 1 n ft 2 ) 
r,fc r, otherwise. 

Note that 

(Gi)- 1 ^) n U{ = { (x',s) I x' eu l k ,se [~sUx r ), t k (x')} }. 



The functions 8 rk are uppersemicontinuous, the proof is similar to the proof of 
Proposition [T5], using the properties of sets satisfying Condition ^. 
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Define the mass transport density a(-). Let a = on R n \(fi 1 U Q 2 ). To define a 
in Q l U Q 2 it is enough to define a on T>i U T> 2 - Let x G T>i\lZi, where Z = 1,2. Let 
y G N"dn t (x) and s = |x — y\. Then we define 



l^ 1 1 



ax 



n-l 



n l (y) 

>S(y) 



V 



where k& are curvatures of dfl l t at Then a satisfies the equation (|20|) on subintervals 
of distance rays of Q l t that lie in X>/, and a = on dT>i and on IZi, I = 1,2. We 
calculate, denoting by K ri j the principal curvatures of dfl l r at the point y r G fij, 
nearest to x 



{— ~ w t )ipdz 

2 Ni 



* J=1 fc=1 -Mn(Gj< )-i(7>0 
1=1 k=l Ju l k 



[ w °G l k -V o G l k ){y o Gi)(*i o G<) J n G< 



x 



2 ATj 



[( s _ r ) - V{G l k {x', 0))] (^(Gi(x', s))J n G l k {x', s) dsdx' 



2=1 k =i M J# r , k v 



) 



da(G l k (x',s)) 
d s 



n-l 



+a(G l k (x', s)) E — V •^(z', s) dsdx' 

2 /V, 

EE L* l k(G l k (x',o)) 
i=i k=i 



X 



2 2V; 

EE / 

l=i fc =i M^)^ 
aDwDcp dz. 



a(G l k (x',s)) d ^ Gk d ^' ,S ^ JG l k {x',s) dsdx' 



a{z)¥ k {z) D(f{z) Dw(z) dz 



R» 

Thus we have proved that the mass balance equation is satisfied. This implies The- 
orem 0. 



VARIATIONAL EVOLUTION PROBLEMS AND NONLOCAL GEOMETRIC MOTION 29 



8. Compression molding model 



Compression molding is the process of deformation of an incompressible plastic 
material between two horizontal plates. The following simplified mathematical model 
of the process was derived by G. Aronsson [ |Ar95| ] based on Hele-Show model for non- 
Newtonian fluid. Suppose that the distance between the horizontal plates is small. 
Then we can assume that the region occupied by plastic at each time t has the form 
Ut — fit x [0, h t ) where fl t C R 2 , and that the pressure in plastic does not depend on 
the vertical coordinate, i.e., pressure is the function u(x,t) where x G fit- Evolution 
of rescaled fit and u(x, t) is described by the following free boundary problem. Given 
an open set f2 G R 2 find an expanding family of open sets Q t G R 2 ,t > 0, and a 
function u(x,t) defined on U t (fi t x {£}) such that 



(60) 
(61) 
(62) 



div(|Du| 



p-2 



V 



Du) = 1 
u = 
\Du\ p - 2 



in fl t , 
on T t , 
on T t , 



where T t = dfl t , V denotes the outer normal velocity of T t . Condition 
that the free boundary T t moves with the velocity of the flow. 

In the paper [|AE|| the asymptotic limit as p —>■ oo in the problem (|60| 



means 



(|62|) was 

considered. This limit corresponds to the case of highly non-Newtonian fluid. 

It was shown in | |AE| | that formally sending to a limit in ( |B"(1D - ( |6"2"| ) one obtains 
the following problem. Find a family {fit} of open subsets of R 2 , a function u(x,t), 
and a mass transport density function a(x,t) satisfying: 



(63) 



\Du\ < 1 
a > 

supp(a) G {\Du\ = 1}, 
— div(a-Dw) = 1 



a.e. in fl t , 
a.e. in fl t , 



in fl t , 

where the last equation is understood in the weak sense, and 
(64) u = 0, V = a on T t , 



where V is the outer normal velocity of T t . 
It was shown in ||AE|| that solutions of (^' 



B3) have the form 



(65) 



u(x,t) = d Qt (x), 



where the right-hand side is defined by ([!]), and that formally the problem 
can be rewritten as following. Find {fit} such that 



(66) 

where 

(67) 



w - d t w G 19/00 [u] 



where x^t(') * s ^ ne indicator function of fl t that equals 1 inside and outside fl t - 
Existence of a weak solution of - ( |5BD was proved in ||AE|| . Namely, there exists a 
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family {fit} of sets of finite perimeter such that d t w is a nonnegative Radon measure 
and 

/ w{x, t)(v(x) — u{x, t))dx < / (v — u{-, t))d(d t w(-, £)) 

(68) 

for every v with < 1, for a.e. t. 



The following law of motion of the free boundary T t was derived in ||AE|| by a formal 
calculation: 

(69) y = 7 (l_^), 

where 

outer normal velocity of T t , 
function 7n 4 (-) defined by (|21~1), 
curvature of I\. 

The equation (^Uj) was derived as following. Starting from fl5"3|), we perform calcu- 



lations similar to the ones that lead from ( |T8"D to (^). Thus we deduce that equation 
(|63|) can be formally rewritten on each distance ray R x , where x G T t , as ODE 



(70) a '( s )_ a ( s )^^ + i = , S 6(0, 7nt (x)). 

1 — 

By the nature of mass transfer process in the compression molding model (i.e., mass 
transfer from within the set onto the boundary), we expect that the mass transport 
density equals to zero on the ridge of Q t and equals to the outer normal velocity at 
the boundary. This translates into the following boundary conditions for the ODE 
dm on R T : 



(71) a(0) = V(x), 0(7(2:)) = 0. 

The function a(-) and the constant V can be found from ([T(]), ([711). V has the 
expression (|69"D. 

In this section we prove the connection between the variational equation (|66| ) and 
the geometric equation (|69"D. 
Let £,T be defined by (§). 

Theorem 30. Le£ {fit}, t G be an expanding, locally Lipschitz continuous fam- 
ily of open bounded sets. Let for every t G the set Q t satisfy Condition [| with 
radius r > 0. Suppose that the equation ffiQ) is satisfied at every point of (2,1) 
differentiability of the surface T. Then the equation fthDi) is satisfied for a.e. t G R , 
where the functions u(x,t), w(x,t) are defined by $6$), W\) - 



Proof. The proof goes along the lines of the proof of Theorem [|. We will sketch 
the proof and present some details for the parts that are different from the proof of 
Theorem [5|. 

The main steps of the proof are following: 
Step 1. Definition and properties of mass transport density. 
Step 2. Show that the main mass balance equation is satisfied. 
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Step 3. Show that the assertion of the Theorem follows from the main mass balance 
equation. 

We discuss each step. 

Step 1. 

Let (x,t) be a point of (2, 1) differentiability of the surface T in R 2 x R 1 . Define 
mass transport density a(y,t) = a(s) on the distance ray R x of fit as the solution 
of flT0|), QTTJ) . Do this for every such point (x, t). It follows from the hypothesis of 
the Theorem and from Proposition |23| that a(y, t) is now defined at C 3 a.e. point of 
the set E = U t (tt t x {t}) C R 2 x R^. Define a(y,t) by zero at all other points of 
R 2 x R^. 

Solving ([TCP, ( [71] ) (and taking into account (j69|) ) we see that on the ray R x 

1 — KS (1 — K7) 2 



(72) a(s) 



2k 2k(1 — ks) 



This function satisfies assertions (a) and (d) of Proposition |26]. Instead of assertion 
(c) of Proposition we have 

(73) \a(s)\ < C(-r(x) -s), \a(s) -V(x)\ <Cs. 

To see this, we note that from the conditions of the Theorem and (^) it follows that 

V{x) < C(r , diam^t). 
Then the assertions (a) and (d) of Proposition [26] for the function a(s) defined by 



([72]) follow from (^) like in the proof of Proposition ^6 
Let us prove (|73|). We can rewrite (^) as 

a{s) = — - — 1 + 



2 V l-KS 

Now, since < s < 7 < diamf2 t and firy < 1 and k > — ^, we get the estimate 
(74) < l^CL < c(r , diamfi t ), 

1 — 

and the first inequality of (|73|) follows. To prove the second, we calculate using (|7T|) 



a(s) —V = a(s) — a(0) = - 



(1~«7) ; 

1 — KS 



and use (^) and the inequality 

< 1 — < C(r , diamfi t ) 

to finish the proof of (|73|). 

The purpose of this step is to prove that for any smooth function ip G C^(R 2 x 
(0, T)) we have 

(75) f [ w(tp + ip t )-aDuD<p = 0. 

Jo Jr. 2 

This is main mass balance equation for compression molding model. 
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The proof follows Section |5]. 

Fix t. We use local coordinate systems on (Vt t ) r where < r < r , defined in 
Section |3] and notation introduced there, in particular sets Uk, U^, partition of unity 
coordinate mappings Gk and functions %^ defined by (p9|), (|30|) , ( p7|) and ( |3T]) 
respectively (k=l,...,N). 



We use notation J5 n T r introduced in fl2T)|). From Proposition |9] it follows that for 
r G [0, r ] 



The function 



r r = {(x,t)GR 2 x(0,T) I dUx) = -r}. 



(x,t) -> d£, t (x 



is Lipschitz since {fit} is a Lipschitz family of sets. Then by Proposition 3.2.15 of 
|Fcd69|| , for a.e. r G [0,r ] 



(76) set T r is (H 2 , 2) recti f table. 

In addition, by Lemma 24, the function — > d,Q t (x) — C\x\ 2 , where C is large 
enough depending on ro, is concave in x for every t in the set -E ro . Then it follows 
from Theorem 1 of Appendix 2 of [Kry_87] and Propositions ^3|(b) and 11 that for 
a.e. r e [0, r ] 

(77) surface T r is (2, 1) dif ferentiable at 7i 2 a.e. point] 

functions dtd^Jx), Dd s nt (x), D 2 d s nt (x) are bounded and 

(78) 7i 2 measurable on T r . 

Function (x,t) — > Jn t ( x ) is uppersemicontinuous since the family of sets {fit} is 



Lipschitz. Then similar to Propositions 27, E9 and Corollary Bq we prove that the 



mass transport density a(y,t) defined at Step 1 is Lebesgue measurable in R 2 x R 1 ^. 
In the proof we use ( |78|) (i.e., we choose such r in the definition of local coordinate 
systems that ([78]) is satisfied). 

It follows from Proposition ^ that for a.e. t we have the following: a.e. point 
x G Qt lies in the relative interior of the distance ray of the set fl t that intersects dQ t 
at a point of (2,1) differentiability of the surface T. Fix such t. Let r G (0,ro]. Let 
y G dQt be a point of (2,1) differentiability of T. Let be the curvature of dVL t at 
let denote be the curvature of d(Q t )r at the unique point y G <9(fi t ) r nearest to x, 
then 



1 1 
— = — + r. 
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We calculate using properties of the function a on the rays proved in the Step 1, 
equation (|70|), and Lemma |T^ with n = 2: 



wip dz 



R 2 



<p dz 



V / (<p o G k ) (% o G k ) J 2 G k dx 

k=l J Uk Jr 



+a[s — r) 



K,. 



1 — K*Sj 



((p(G k (x r , s))J 2 G k (x', s)dsdx' . 



In the last expression we integrate by parts and use Q56|) , ( |71| ) and (|58|). Then we get 



[ wipdz = V / * fc (G? fc (a/,r)) 



fc=i 

*r*W ^ d(p(Gk{x' , s)) , 
+ / a(s-r) J 2 Gk{x ,s)ds 



a{0)ip{G k {x',r))J 2 G k {x',r) 

dx' . 



ds 



V (Pt(y))d t d s nt (P t r (y))(l-r Kr (y))dH 1 (y) + / aD^Dudz, 



where 
(79) 
and 
(80) 



d(n t ) r = {y E d(n t ) r > r], 



Pi : (a)r -> fit 

is the nearest point projection mapping. By Proposition |9] b) the mapping P r * is 
well-defined and onto. 

By Proposition [3~3|, if T r is (2,1) differentiable at y G d(fi t ) r , then T is (2,1) differ- 
entiable at y' = P*(y) and 

dtd^PXy)) = d t d s nt (y) 

By the equality 

d fa t ) r ( x ) = d k( x ) + r for (M) e ^ > >" e (0,r ] 

we see that 

^(^(y)) = ^ )p (i/) 

for y e 9(fit) r . 

Thus we have showed that for a.e. t € [0, T], every r e (0, r ] the equality holds 

(81) / wtpdz = I {^oP^dtdl^il-rnDdH 1 + / aDtpDudz. 
Jn 2 Jd(n t ) r J &t 
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In the next lemma we show that in the first integral at the right-hand side of (|8^ 
we can integrate over the whole boundary d(Q t )r- 



Lemma 31. Let r G (0,r ] be such that properties ff7b\), (|7%), ( |7^ are satisfied. 
Then for a.e. t G (0,T) 



Id( 



ipdtd^il-rKDdH^O 

'd(n t ) r n{-y (nt)r <r} 
Proof. We will show that for a.e. t 

(82) d t d s [nt)r (y)(l - rKt{y)) = H 1 a.e. on d(Q t ) r C {j { n t)r < r}. 

By ( f77|) for a.e. t the surface T r is (2,1) differentiable 7Y 1 a.e. on <9(f2 t ) r . Fix such 
t. Then in ([82] ) we can consider only y at which T r is (2,1) differentiable. Fix such y. 
By Proposition |9|, 

7(n t )r > r on d(tt t ) r . 

Thus 

7(a*)r(2/) = r - 

It follows that 

<{y) < -• 

r 

Let y' = Pr(y), i-e., y' is the point of 3VL t nearest to y. Consider 3 cases. 

Case 1. y is the unique point of d(Q t )r for which y' the nearest point on dVL t - 
Then the calculations of [ EGan , Proposition 7.1, Steps 2,3 imply that 

4(v) = -■ 

r 

Thus we have (|82|) in this case. 

Case 2. There exists z G d(Qt)r, z ^ y, such that y' = F£(z), and 

B r (y)nB r (z)^$. 

Note that 



\y-y\ = 

Then, denoting w the point ^r, we get 



\ z — y 



We have 



w G B r (y) U B r (z). 



B r {y), B r (z) cR 2 \n t . 



Thus w G R 2 \ fit, and y' is the point of dVL t nearest to w. Denote 



v = y' + 



\w — y 

Then, by Condition [| we get the following: 

v ^ Q t , y' is the unique point of dfl t nearest to v. 

But 
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Thus we proved that 

v G <9(f2f) r . 

The points y and z divide the circle dB r (y') on two arcs, and the point v is the 
middle point of one of these arcs. Denote this arc C\. Repeating the above argument 
inductively and using continuity of distance function, we prove that 

Ci c <9(O t ) r . 

But then 

4(v) = -■ 

r 

Thus we have ( ]8"2]) in the Case 2. 

Case 3. There exists z G d(Q t ) r , z ^ y, such that y' = Pf(z), and 

B r (y) n B r (z) = 0. 

Then 

dB r (y)ndB r (z)=y'. 
Since {fit} is an expanding family of sets, we have 

(83) B r (y), B r (z) cR 2 \fl r for r < t. 
If 

(84) y' G dQ T » for some r* < t, 
then the same is true for all r G [t*,£], and ([8*3]) implies that 

y G d(Q T ) r for all rG[r*,t]. 

Then 

^) P (l/) = 0. 

Thus we have (|8~~"f) if (|84"|) is satisfied. 
The remaining case is 

(85) y' £ dtt T for all r < t. 
Introduce a coordinate system (x^a^) on i? 2 in which 

y' =(0,0), y=(0,r). 

Then 

z = (0, -r). 

Since {fi t } is a continuous expanding family of sets, the function 

0(r) =dist(y , ,fi T ) 
is continuous and nonincreasing, and 

0(t) = 0. 

Let w T eAfanM- Then by © 

(86) w r G 0B*( T) (O, 0) \ [S r (0,r) U B r (0, -r))] . 



36 MIKHAIL FELDMAN 

In particular, 

(87) w T — > y' as r — > £. 
By Condition |3] we get 

(88) -S r (v T ) CR 2 \H T where v T = w T + Ay' - w T ). 

\y'-w T \ 

By (p6|), ( p7|) we see that there exists a sequence r,- — > £ such that 

f Tj — > p, where p is either (— r, 0) or (r, 0). 
Then by (|88|) and continuity of the family {Qt} we conclude 

either B r (-r, 0) C R 2 \ H t or B r (r, 0) C R 2 \ tt t . 

Let in fact 

B r (r,0) cR 2 \n t . 
Then y' = (0, 0) is the point of dVL t nearest to the point (r, 0). Thus 

(r,0) E d(Q t ) r . 

But y = (0,r), and thus 

B r (y)nB r (r,0)^d). 
Thus the points y and (r, 0) satisfy the conditions of Case 2. Thus we get 

<(y) = -• 

r 



Case 3 is proved. Then Lemma [31] is proved. □ 

The properties fl76|), ([77]), © are satisfied for a.e. r e (0,r ]. Then there exists 
a decreasing sequence — > 0, where i = 1,2,... and G (0, r ], such that (|76l), 
(|78|) are satisfied for each ty Note that by (|78"D , the function 

(y,t) -> 4.(y) 

defined on r r . is 7-^ 2 measurable. 



Now we can integrate ([81]) by t and use Lemma |3l] to get 

[ T [ w V dzdt = [ T [ (voptJddtn^Xl-^lJdH'dt 
Jo Jn 2 Jo Ja{Q t ) ri 

(89) + f [ aDtpDudzdt. 

Jo Jn t 

It follows that the first integral at the right-hand side does not depend on ?y Thus 
it is enough to compute the limit as i — > oo. We will prove the following: 

(90) lim / / (^0^)94^1-^)^1=- / / XnA^dxdt 

Jo Jd(n t )r t Jo Jn? 

We first prove that such equality is true if the boundary satisfies additional regu- 
larity assumptions. 
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Lemma 32. Let {f^} satisfy fffflj) - fffflj) with r = 0. Let V(x,t) be the outer normal 
velocity ofT t , defined by (|6[) at every point (x,t) of differentiability ofT and defined 
as at all points where V is not differentiable. Then for every ip G C*(R 2 x (0,T)) 

(91) f [ tpVdH 1 dt = -[ [ xnA^dxdt 

JO JdQt Jo Jr. 2 



Proof. Let the set E be defined by @. Then by (|7q) E has locally finite perimeter. 
By© 

H 2 (dE\d*E) = 0, 

where d*E is the reduced boundary of E. Let $ : R 2 xR 1 ^ R 2 x R 1 be defined by 

$(x,t) = (0,...,0,<f(x,t)) 

Let i>E(x,t) be the measure-theoretical outer normal to T at (x,t) G T. Then by 
Green-Gauss theorem for sets with finite perimeter ( |[EGar92|l , section 5.8) 

rT p 

Xn t dt¥>dxdt= I div$ dxdt = I ^UEdTi 1 



o Jr 2 Je Jr 

where 

div$ = — — -, where X3 = t. 

-f-f OXi 

1=1 

At every point (x, t) of differentiability of T we have 

u E (x,t) = -D {X:t) d s E (x,t), 

where d s E is the signed distance to the boundary of the set E in the (x, t)-space. Thus 
we get 

! d t <pdxdt = - [ ipd t d E dH 2 . 
Je Jr 

Let / : r — > R 1 be the mapping defined by 

f(x,t)=t 

Then / is Lipschitz, and / and T are differentiable at Ti 2 a.e. point (a;, t) G V. At 
such point (x,t) the gradient of / is a linear mapping Df(x,t) : T X}t T — > R 1 , where 
T x>t T is the tangent to T at (x,t) space. Let e±,e2 be such orthonormal basis in R 2 
that C2 is the inner normal to dVt t at x. Let r be the unit vector in the t-direction. 
Then the vectors 

ei, e 2 = -d t d E (x,t)e 2 + \D x d s E (x,t)\r 
form an orthonormal basis in T x t T. We calculate: 

Df(x,t) ei = 0, 
Df(x,t)e 2 = \D x d s E (x,t)\. 

Thus 

\Df(x,t)\ = \D x d s E (x,t)\. 
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Also, the following relation holds: 

D xt d s ixA) = — — —D xt d s F (x,t). 

' QtK ' \D x d" E (x,t)\ ' Ey 1 

Now, applying formula 3.2.22 of [ Fed69|l (which is applicable by (|76|) ), we get: 
J^d t d s E dH 2 = J^d t d s Qi \Df(x,t)\dH 2 

[ ^d t d s nt dn l dt 
I LpVdH l dt. 

Jdn t 

The lemma is proved. □ 
Now we can prove (|90|). Each {{Q t ) r .} satisfies (|76| ) - ([781). Thus we have 

/ / V V^dH 1 dt = - [ [ X (n t)r dtVdxdt, 
Jo Ja(n t ) ri Jo Jr 2 



(92) 



where V Tl is velocity of (f^) rr It follows from Proposition that C 3 (T) = 0, and 
thus 

X(Q t )rS x ) Xn t (x) for a.e. (x,t) G R n x (0,T). 

Then by Dominated Convergence Theorem the right-hand side of (^) converges to 
the right-hand side of ( |90| ) as i — > oo. Thus it remains to prove that the left-hand 
side of ( |92| ) and the left-hand side of ( |90| ) converge to the same limit. 

Let P* r . : d(fl t )r — * d(Qt) ri be the nearest point projection (well-defined by Propo- 
sition |9|). P* r . is a Lipschitz map by Propositions^, [II]. Then using Lemma H and 
identity (p3|), the difference between the left-hand sides of (^) and ( p0| ) can be trans- 
formed to 

(93) f [ (ipo Pl n - <p o Pi) V ri o Pl n (1 - TK^dH 1 dt + 

Jo Jd(n t ) r 

n f T [ <p ° ^ ° 4 dn l dt = i u + Ti i u . 

Jo Jd(n t ) r 

We have \V n \ < C, where C does not depend on i, and < -. Using the fact 
that P r * r . — * P r * as fj — > and Dominated Convergence Theorem we see that Ii ti — ► 0. 
We also have |/ 2i j| < C. Thus the expression ( p3|) converges to zero. Thus (0) is 
proved. 

The equalities fl8g) and imply fl7§). 

Step 5. 

The family {fit} is expanding. Thus the left-hand side of ( |9"U] ) defines a nonnegative 
linear functional of <p G C°°(R n x [0,T]). By ||LCar92|| , Chapter 1.8, Corollary 1, it 
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follows from fl90| ) that d t Xn t 1S a nonnegative Radon measure. Thus the mass balance 
equation (|75| ) can be rewritten as 

/ wip dxtdt + I cp dwt — / / aDuDip dxtdt = 0. 

'r 2 Jr 2 x(o,t) Jo Jn 2 

From this equation it follows that for a.e. t G [0, T], every ip G C*(R 2 ) 

w(- 1 t)(fdx+ / (pdw t (-,t)— / a(-,t)Du(-,t)D(p dx = 0. 

R 2 JR2 JR2 



Now, repeating argument of Section y, we conclude the proof of Theorem |30|. □ 

Appendix A. Appendix 

Let fl be an open set. We write d(x) for d^x) below. The purpose of this section is 
to prove that the gradient of d^ is locally Lipschitz at x G fl \ 1Z and give an estimate 
of the Lipschitz constant in the terms of the distance between x and endpoints of the 
distance ray R x . 

Proposition ^ should be compared with Proposition 4.1 of ||EGan|| and with in- 
equality 4.8(8) of |edS|. 



In the inequality (^) below the only assumption regarding the point X\ is that it 
is close enough to x. In particular it is possible that x\ G 1Z. In the inequality 4.8(8) 
of | |Fed59j j the conditions on x and x\ are symmetric and exclude the possibility that 
x x G K. 

Proposition |33] improves the estimates of Proposition 4.1 of ||EGan|| in the following. 
Two quantities are estimated explicitly in Proposition [J3[ the local Lipschitz constant 
of Dd(-) at x and the size of the neighborhood of x in which the estimates (|95|) and 

am hold. 



Proposition 33. There exist constants C and M depending only on n such that the 
following is true. Let Q G R n be an open set. Let e > 0, x G Q\7Z, and let 

(94) d n (x) > Me, 7n (z) - d n (x) > Me. 

Then by Remark^ there exists a unique ray R x . Denote y and v the lower and upper 
ends of R x , i.e., y = R x fl dQ and v — R x fl 71. Then for every x\ G VL satisfying 
| x — x\ \< e, the inequality holds 

(95) \y-yi\<c(l+ \ x ~ y ^ \ \ x-x 1 \, 



, x — v i 

where y\ G Mm{x\). If in addition the function d(-) is differentiate at x\ then 

C 

(96) \Dd(x) - Dd{xi)\ < —\x - x x \. 

Proof. We first prove the inequalities (p5]), (p^) assuming that 

(97) \ x — v \>\ x — y \ . 

Denote d := du(x). By fl9"T|), we can find a point O on the interval of R x between 
v and x such that \ x — O \= d. Let y G Afan(x), then we also have y G J\fan(0). 
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In the calculations below C will denote different constants depending only on n. We 
assume that 

(98) M > 10. 

Choose X\ G B e {x). Let X\ be projection of X\ onto R x , then 

I ii - - 1Q 

by fl94|), (|9~8|). Thus x\ lies between y and O on R x . Let di = \x± — 0\, then 

(99) r ^i^w 

Introduce the following coordinate system in R n : let the point O be the origin, let 
e n — i^Iqi (thus e n is the unit vector along the ray R x ), let e%, e n _i, e n be an 
orthonormal basis in R n . Then in these coordinates 

z = (0,d), fi = (0,di), y = (0,2d), 

where G R n_1 . We also have c?i), where x' x G R n 1 . Let ?/i G <9f2 be 

such that |xi — yi\ = d(xi). Let the coordinates of y\ be (yi,yi, n ) where y[ G R n_1 , 
yx, n e R 1 - Then 

(100) \x x - y x \ < \x x - y\ 
or 

(101) 1^ - y[\ 2 + |di - yi , n \ 2 < \x[\ 2 + \d x - 2d\ 2 

Claim 34. There exists M 1; depending on n, such that for any 5 G (0, 2) the follow- 
ing is true: 

Mi 

(102) if d > — j-e then \y[\ < 5d. 

cr 

Proof. Note that 

(103) K\<e, \y[\<2d. 

Indeed, the first inequality is true since x\ G B £ (x). To prove the second inequality 
of flTUD we use (HOI) and (II) to S et 



|yi - x[\ 2 < e 2 + \d + e\ 2 < (d + V2e] 

so 

(104) \y[\ < \x[\ +d+ V2e < d + 3e < 2d 

since d > lOe. Thus ( pL03| ) is proved. 

Suppose that the assertion of Claim |34| is false, that is 

(105) \y[\>5d. 
By (HD 

Mi 

(106) 5d > — -£ > 2e 

o 



2 
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if Mi > 4. Then we get from ([101]), QID5D , @, QIOSD 



(Z/l,n-rfl) 2 < £ 2 + (2rf- C / 1 ) 2 -(K|-|x / 1 |) 2 

(107) < £ 2 + {^d) 2 -{5d-e) 2 

S 2 

< 4d 2 (l )-2d 2 + 25de 

< 4rf 2 (l-^) 

The last inequality follows from —2c? 2 + 25de < which holds by ( |102| ) if 5 < 2 and 
Mi > 8. 

Consider two cases: 

Case 1 : y l n > 6?!. Denote by cv the angle between vectors yi — X\ and (x' 1 ,yi tn ) —x\. 
In the triangle X\,y\, (x[, y\, n ) the side (x[, y\, n ) — Hi is orthogonal to X\ — (x' 1 ,yi jTl ), 
since (x[, y ljn ) - yi = (x[ -y[,0), and x x - {x[, yi, n ) = (0, d 1 - y hn ). Thus we get 



tana 



m ~ x i I 
yi, n - d\ 



But 



Iz/i — x '\\ > Hl/il _ kill = \8d — e\ = 5d — e 
by (TO . By ([lO^) we get 



tana > 



5d — e 



2djl -f 2*/l , 



= (l--" £ )> 



5 



4a/1 



6 2 



where the last inequality follows from ( 106 ). Thus, using the condition 5 G (0,2) we 
get 



(108) 



cos a 



Vl + tan 2 



< /4(4- i! ) < 
16 -35 2 " 



2(16-35 



2^ ' 



Consider the triangle xi,yi, (x[,0). The angle between the vectors (x^O) — x\ 
and y\ — x\ is 7r — ft, and so we get 



(109) 



\y x - (xi,0)| 2 = \x\ - K,0)| 2 + \y x -£i| 2 - 2\x x - (xi,0)||j/i - x x \ cos(vr - ft). 



Since x\ = (x[,di), we have 
From ( |101| ) we get 



\Xi — [x 



ij 



d\. 



\yi - xi\ < v^ 2 + (2d - rfi) 2 . 
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Then by ([T09f ), ( |T08|) we get 

bi-K,0)| 2 < dj + [e 2 + (2d - di) 2 } + 2d iv /e 2 + (2d - dif (1 - ^) 

< d l + (2d _ + 2di(2d - dO(l + 2(2 /_ di)2 )(l - 33) + ^ 

< 4d 2 + d 17 - ^— + e 2 -^-d 1 (2d-d 1 )6 2 . 

2d — di 16 

Thus we see that 

(110) |y 1 -(x' 1 ,0)| 2 <(2d-e) 2 , 
provided 

(in) d nTT—r - ~k d ^ 2d - d ^ 2 < ~ Ad£ - 

2d — di 16 
Using ( |99| ) we see that ( |111| ) is satisfied if 
(112) d 2 5 2 > C(e 2 + de) 

where C is a large enough constant. In its turn, ( |112| ) is true if d > ^rr£, where Mi 
is large enough. Thus, for such d, the inequality ( |110| ) holds, and so 

\y 1 -(x' l ,0)\<2d-e. 

Then 

\yi\<\x' 1 \ + \yi-(x' 1 ,0)\<2d=\y\. 

This contradicts the fact that y £ A/an(0). Thus ( |105|) is false. Thus Case 1 of Claim 
3~4| is proved. 

Case 2: y l n < d%. 

Using this condition, we get 

\(x[,0) - yi | 2 = |(xi,0) -xi +xi — 2/i| 2 = I - (0,di) + (x[ - y[, d 1 - yi, n )\ 2 
= dl + \xi - yi\ 2 + 2(y hn - di)di 

< d\ + \xi-yi\ 2 <d\ + e 2 + (2d-di) 2 
= Ad 2 - 2di(2d- di) + e 2 < (2d-e) 2 , 

if di(2d-di) > 2de. 

The last inequality is in fact true by and d > lOe. Then, as in the Case 1, we 
arrive to a contradiction. Thus Case 2 cannot happen for d and e defined in Claim 
34. Thus, Claim 34 is proved. □ 

Now the rest of the proof of part (a) of Proposition [IT] follows the proof of Propo- 
sition 4.1 in IIECan! . 
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The main step is to show that there exist constants C and M depending only on 
n such that the inequality d > Me implies 

(113) \y - 2/x | < C\x 1 - xt\. 

Let S G (0, 1) be chosen later, and let d > ^-e where M\ is the constant from 
Claim |34j. Then by Claim we get 

(114) \y[\ < 5d. 
Since y G N"au(0), we get \yi\ > \y\ — 2d, and so 

\y[\ 2 + \ yi , n \ 2 >4d 2 . 

Thus we estimate using ( |114|) 



(us) > 2d -Ml!- C Mr> M _MJ! (1+cn 

where all constants C do not depend on y[ and d. So 

12 



< 2d - d x = y hn -dt + 2d- y 1>n < {y 1>n - d x ) + ^-(1 + C5 2 ). 
Now we can estimate 

(2d-d 1 ) 2 < (y 1>n -d 1 ) 2 + (y 1>n -d 1 )^(l + C5 2 )+^(l + CS 2 ) 

(116) < (y 1 , n -d 1 ) 2 + (y 1>n -d 1 )^(l + C5 2 )+C5 2 \y' 1 \ 2 

where we again used ( |114| ) . Using ( |101| ) , (^) and ( |103| ) we get the following estimate 

e 2 



Vi,n-dx\ < ^e 2 + (2d-d 1 ) 2 = (2d-d 1 ) x l 



(2d-d 1 ) 2 



e 2 , .... . „£ 2 . 



< (2d-d 1 )(l + C5 4 ) since d>^-e. 

o 2 

Thus, by ( TO ): 



(2d-d 1 ) 2 < (y hn -d l ) 2 + (2d-d l )^-(l + C5 4 )+C\y' 1 \ 2 5 2 
(117) < (y 1 , n -d 1 ) 2 + Q\y' 1 \ 2 , 



where 



Using (|99| ) we estimate 



= (1-^)(1 + C5 4 ) + C5 2 . 



e<^(l + C5 4 ) + C6 2 <l, 
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if 5 > is small enough. Fix such 5. Let M\ be the constant from Claim |33| defined 
by this 5. Let M = ^ be our choice of the constant M in the inequalities (p4|). 
Then the inequalities ( |9lD imply ( 117 ) with < 1. Thus ( |117|) and ( |101| ) give: 

(118) \y'x-x'i?<K? + ®\y'x?- 

Now we estimate: 

1 2/i |<| y[ -x' x I + I x' x I 

and from this, using ( |H8|) , we get 

I 2/1 T < |^-x' 1 | 2 +2|^-a; / 1 ||a; , 1 | + |x / 1 | 2 

< (i + / i)(Ki 2 + eKi 2 ) + (i + -)|x;i 2 . 

n 

Choose fi > so small that (1 + //)©< 1, then we get 

(119) \y[ |<C|xi \=C\x l -x 1 | . 

It remains to estimate | yi iTl — y n | where y n = 2d. 
From ( |115| ) we get 

(120) 2 d-y 1>n <C^-<C\y , 1 \, 
since | y[ |< 5d. 

Since x± = (x\, di), d(xi) —\ X\ — y\ |, and d(xi)) —\ x\—y \= 2d — di, we calculate 
(using the fact that Lip[d(-)] = 1) 

2/i,n — di < \ yi — x\ \= 2d — di + d(x\) — d(x\) < 2d — d\+ \x\ — X\\ 

and so 

(121) y hn - 2d <\ xi - xx \ . 

The estimates O, CIT20D , and O imply ( |TT3l) . 
Now since 

^ — 2/i n ,/~ x ^1-1/ 
Ddixi) = -. r , Ddixi) = — t, 

\xi-yi\ \x 1 —y\ 

we get the following estimate 

\Dd{x!) -Dd(x x )\ < 

then using the equality \xi — y\ = 2d — d\ and the inequalities (^) and (|117| ) we 
deduce that 

C 

\Dd(x\) — Dd(xi)\ < ~~^\ x i ~ 

Note that Dd(x\) = Dd(x) since the points x and Xi are on the same ray R l x . Now 
from \x\ — x\\ < \x — Xi\, we get (^) in the case of (P7|). Note also that ( |113| ) and 
inequality \x± — x±\ < \x — X\\ imply 



\Xi — X\ 


+ 


2/1-2/1 




Xy 


-y 
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In the case | x — v \<\ x — y | we can make the following reduction to the case fl9?p. 



Let 
(122) 



Q 



x 



y | — | x — v 



and let fl be the set {z £ Q | dist(,z, dQ) > Q}. Then B £ (x) C Cl since the function 
dist(-,<9f2) is Lipschitz with constant 1. For z G Cl we have 

(123) dist(z, dO)) = dist(z, d&) + Q. 

To see this consider z\ G J\fdn{z), and let z 2 be the point of intersection of the interval 
connecting z and Z\ with dQ. Then z\ G Afaui^) since Z2 lies in the relative interior 
of the distance ray R Zl that intersects dfl in z\. Thus | Z\ — z 2 \— Q. Now suppose 
that there exists a point z 2 G dQ such that | z — z 2 \<\ z — z 2 |. Let z\ G Mq^.[z 2 \ 



then \z\ — z 2 



|= Q. Then 

z — Z\ I < I z 



Z2 



-Q 



Z — Z\ 



\ z 2 - Zi \<\ z- z 2 

a contradiction to the fact that z\ G Mao, (z). This proves (|123|) . So the distance rays 
and the ridge set of the set Q are intersections of those of the set Q with the set Q. 

Then in ( ^B| ) we can consider d(-) as the distance to dfl, so we can consider Q 
instead of Q, and then the inequality ( p7|) is satisfied. 

To show (^), denote by w and Wi the points of intersection of dCl with the intervals 
(x, y) and (xi,yi) respectively Then it follows from (|123|) if G A/" 9 q(x) and -tfi G 
and 



w 



y 



wi - y% |= Q. 



Then it follows from ( |122| ) that 



x — v 



x — w 



Using ( |122|) and letting A 

y = x + A{w 



x—y\ , 

— 4, we get 



Xi + 





— If 




- 1 



(W! - Xi). 



The inequality fl95|) can be applied to the set f2, so 



10 



w l \< C 



Now we estimate: 

I y-yi 



< 
< 



X 

+ 

C(l 
C(l 



Xi + A 



X 



X 



Xi 



x 



w 



X 



w 



X\ — Wi 

hA)(\x- 



A) 



X 



I Xi — W i 
- 1^ (^1 - Xi) 

• Xi I + I 10 — Wi I) 
Xi I . 



-{wi - Xi) 



This implies 



□ 
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